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A method for deriving quantum kinetic equations with initial correlations is developed on 
the basis of the nonequilibrium Green's function formalism. The method is applicable to a wide 
range of correlated initial states described by nonequilibrium statistical thermodynamics. Initial 
correlations and the real-time evolution are treated by a unified technique employing many- 
component "mixed" Green's functions. The Dyson equation for the mixed Green's function 
leads to a set of equations for real-time Green's functions and new (cross) components linking 
initial correlations with dynamical processes. These equations are used to formulate a generalized 
Kadanoff-Baym ansatz for correlated initial states. A non-Mar kovian short-time kinetic equation 
is derived within the T-matrix approximation for the self-energies. The properties of the memory 
kernels in this equation are considered in detail in Born approximation for the T-matrices. 
The kinetic equation is demonstrated to conserve the total energy of the system. An explicit 
expression for the time-dependent correlation energy is obtained. 



I. Introduction 

The problem of memory effects and initial correlations in nonequilibrium systems is as 
old as transport theory. In 1960s, some general aspects of this problem were studied on 
the basis of fundamental principles of statistical mechanics, and formally exact equations 
for the one-particle distribution function have been derived, involving initial correlations 
and memory HI - 0. However, since a many-particle system "forgets" irrelevant details of 
its initial state, the short-time memory effects and the evolution of initial correlations were 
long thought to be purely theoretical problems related to the justification of Boltzmann- 
like kinetic equations and hydrodynamic equations. Therefore, no serious attempts were 
made to construct explicit kinetic equations describing the early stage of evolution of a 
nonequilibrium system. Now this topic has come to more practical importance due to 
experimental studies of fast relaxation processes caused by the interaction of very short 
laser pulses with matter |^. Another field of high interest in short-time quantum 
kinetics with initial correlations are nuclear collisions 1511 - HI. 
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Although some interesting facts showing the interplay between collisions and corre- 
lations in a many-body system were noticed many years ago (see, e.g, [|10|), a detailed 
investigation of short-time kinetics with initial correlations has been undertaken only re- 
cently by Kremp et al. |]Tl], using the density matrix approach. In these works, the 
quantum BBGKY hierarchy for the reduced density matrices is truncated at the three- 
particle level. This allows one to consider the dynamics of two-particle correlations on the 
same footing as kinetic processes described by the one-particle distribution. Numerical 
calculations |]T2[ show that if the truncation procedure is consistent with some necessary 
conditions, say the conservation of the total energy of the system, then the density matrix 
approach gives physically reasonable results for the time behavior of the Wigner function 
and the correlation energy. The advantage of this method is that one is dealing with the 
single-time reduced density matrices which are the natural ingredients of kinetic theory. 
On the other hand, using the truncation procedure, one has from the beginning to work 
with approximate equations. Then the problem arises how to justify and improve these 
equations in a systematic way. 

It is well known that the real-time Green's function formalism, which to a large extend 
is based on the fundamental ideas of Kadanoff and Baym and Keldysh provides 
an alternative to the density matrix approach to kinetic theory. Within this formalism, 
the kinetic equation can be derived as a component of the exact Dyson equation on the 
Keldysh directed contour in the complex time plane. This is one of the advantages of 
the Green's function formalism since the structure of the kinetic equation follows directly 
from the Dyson equation and all approximations are only introduced for the self-energy 
which may be regarded as a result of partial summation of the BBGKY hierarchy. In 
the literature several generalizations of the Kadanoff-Baym-Keldysh formalism to arbi- 



trary initial states are available. The most successful attempts are due to Hall pS] , \L6 



Kukharenko and Tikhodeev [|r^, Danielewicz Wagner and Semkat et al. pH 
Physically, all these approaches are equivalent, but they differ in the description of initial 
states and in mathematical formulation. Up to now the main goal of such theories was 
to derive generalized Kadanoff-Baym equations or to construct a modified Keldysh dia- 
gram technique for the Green's functions in the presence of initial correlations. However, 
transport theory is conveniently formulated in terms of a kinetic equation for the Wigner 
function. The question then arises how to go over from the system of equations for the 
double-time Green's functions to a closed kinetic equation for the single-time Wigner 
function. This question arises in any version of the real-time Green's function formalism 
and is known as the reconstruction problem or the problem of ansatz. The task is to 
express the double-time correlation functions in terms of the Wigner function. Within 
the framework of the standard Kadanoff-Baym formalism (i.e., for non-correlated initial 
states), Lipavsky et al. |^ were able to derive exact integral equations which allow one, 
in principle, to solve the reconstruction problem by iteration. Based on these equations, 
they also formulated a simple, but rather successful relation between the time correlation 
functions and the Wigner function, which is called the generalized Kadanoff-Baym ( GKB) 
ansatz and now is extensively used in quantum kinetic theory (see, e.g. [^). Clearly, for 
applications of the Green's function formalism to short-time quantum kinetics the recon- 
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struction problem should be solved with taking account of initial correlations. This is, 



however, a difficult task since the derivation of the GKB ansatz I^TJ rests directly on 
the Dyson equation for the path-ordered Green's function, which in turn implies validity 
of Wick's theorem. For arbitrary initial states. Wick's theorem is not valid so that the 
real-time Green's function does not obey the Dyson equation. The possibility of deriving 
the Dyson equation for arbitrary initial states by introducing path-ordered Green's func- 
tions on an extended Keldysh contour was noted by Danielewicz [0] and then used by 
Wagner to formulate a diagram expansion of such Green's functions. Unfortunately, 
this scheme has not been worked out to an extent that explicit kinetic equations can be 
derived. 

The purpose of this paper is to develop a Green's function approach to short-time 
quantum kinetics, which applies to arbitrary initial states and provides the basis for the 
derivation of explicit kinetic equations. We shall show the theory "in action" and give 
some examples of such kinetic equations. The important ingredient of our approach is 
the extension of the GKB ansatz to correlated initial states. 



The paper is outlined as follows. In Sec. II, we touch briefly on the description of initial 
correlated states and discuss a link between this problem and nonequilibrium statistical 
thermodynamics. We also introduce the thermodynamic Green's functions describing 
initial correlations. 

In Sec. Ill, the path-ordered mixed Green's functions are defined on the deformed 
Keldysh contour in the (x, t)-plane, where the x- variable is associated with the "imaginary 
evolution" caused by initial correlations and the t-variable corresponds to the real-time 
evolution governed by the Hamiltonian. The mixed Green's functions have different com- 
ponents depending on the branch of the contour on which the arguments are situated. In 
addition to the thermodynamic Green's functions and the real-time Green's functions, the 
formalism involves new auxiliary objects, namely the cross Green's functions with differ- 
ent types of evolution. The introduction of these functions is necessary to ensure that the 
full mixed Green's function obeys the Dyson equation on the extended contour. In many 
features our scheme is similar to those considered by Danielewicz |]18[ and Wagner |T9 



However, we do not refer directly to the diagram technique because the Feynman rules 
for evaluating Green's functions depend on the initial statistical ensemble describing the 
system. This is why we prefer to formulate the theory in an algebraic language, where 
the self-energies are related to many-particle Green's functions. 

In Sec. IV the matrix Dyson equation for the mixed Green's function is used to derive 
a coupled system of equations for the real-time and cross Green's functions. In particu- 
lar, we obtain generalized Kadanoff-Baym equations including contributions from initial 
correlations. These equations are shown to be equivalent to analogous equations derived 
previously by diagram expansions ||18| . 

Section V is concerned with the problem of reconstruction of two-time correlation 
functions from the Wigner function. In the case under consideration, this problem is 
closely allied to the evolution of initial correlations. Starting from the Dyson equation 
for the mixed Green's function, we derive a set of exact integral equations providing a 
generalization of the reconstruction procedure developed by Lipavsky et al. . 
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In Sec. VI we consider the T-matrix approximation for the mixed Green's functions 
and obtain explicit expressions for the components of the self-energy in terms of the T- 
matrices. An essential point is that there exist exact relations between the cross compo- 
nents of the T-matrix and its real-time components. These relations are used to eliminate 
the cross components and derive a generalized optical theorem including the contribution 
from initial correlations. 

Then, in Sec. VII we derive for the one-particle distribution function a short-time 
kinetic equation in the T-matrix approximation. The collision integral shows retardation 
and contains extra terms due to initial correlations. As a special case, the properties of the 
collision integral and the spectral function are investigated in Born approximation for the 
real-time T-matrices. We show that the kinetic equation conserves the total energy and 
find an expression for the time-dependent correlation energy. Our results are compared 
with results obtained within the framework of the density matrix method. 

II. Description of Correlated Initial States 

To put our analysis into more straightforward language, we shall consider a non- 
relativistic quantum system with the Hamiltonian (we are using units in which ^ = 1) 

H = + H' = ^ I dr^ ViV^t(ri) ■ V^r^) 
2m J 

I r 

+ - y (iri dr^ dr[ dr'^ {r'y^lVlr^r^) ^p\r'2)^p\r[)ip{r^)ip{r2), (2.1) 

where r- denotes a complete set of single-particle quantum numbers, say the position 
vector Tj and the spin variable a^. Integration over r implies summation over discrete 
quantum numbers. The field operators satisfy the usual commutation relations 

'ilj(r)'ijj\r') — ri'ijj\r')ip(r) = 5{r — r'), ^p(r)'ijj(r') — ri^p(r')ip(r) = 0, (2-2) 

where the parameter rj is equal to —1 for fermions and for bosons. From here on the 
delta function 6{r — r') includes the Kronecker delta with respect to discrete quantum 
numbers. If, for instance, r = (r, a), then 6{r — r') = 5(y — r') 5^^,. 

The second term in the Hamiltonian ( |2.1|) describes pairwise interactions between the 
particles. For simplicity, we shall assume that the potential f (Ir^^ — r2|) does not depend 
on spin. In that case the interaction amplitude has the form 

{r'AlVlr^r^) = v{\r^ - r^l) 6{r^ - r[)6{r2 - r'^). (2.3) 

In general, the Hamiltonian contains additional terms describing interaction of the system 
with external fields. Taking the Hamiltonian in the form (|2.1|) , we thus restrict ourselves 
to situations where the system is prepared in some initial (nonequilibrium) state at time 
and we are interested in the evolution of the system at times t > t^. However, the 
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theory can easily be generalized to the case where the system is subjected into a strong 
alternating external field. 

In formulating the short-time kinetics, the first thing one must do is to specify the ini- 
tial state of the system. This can be done most directly by a full many-particle statistical 
operator ^(tg), where is an initial time. Another possible way to describe the initial 
state is by the reduced n-particle density matrices f?i2...n(^o)- r-representation they 

are defined as 

(^1 ■ --rnlQu-MK ■■■<)= Tr {^\r'J ■ ■ ■ ^\r[)^P{r,) ■ ■ -^(rj g{to)} . (2.4) 

Mathematically, these descriptions are equivalent if all reduced density matrices are given 
and both of them are used in the Green's function formalism with initial correlations. It 



is interesting to note that, in their pioneering work |]T3|, Kadanoff and Baym start with 



the correlated state described by the equilibrium statistical operator 

g{t,) = g^^ = e-Pif^-^^) /Tre-^(^-^^), (2.5) 

where j3 = 1/T is the inverse temperature and fi is the chemical potential. Here the 
correlation effects are incorporated through the interaction term in the Hamiltonian H. 
However, Kadanoff and Baym did not consider the evolution of correlations since their task 
was quite opposite. To eliminate the influence of initial correlations, they take the limit 
— s> — oo and assume that the nonequilibrium evolution of the system is caused by some 
auxiliary external fleld acting on the particles. Further the most part of works devoted to 
quantum dynamics with initial correlations, except for the works of Danielewicz [jl8| and 



Wagner ||T9[ , was based on the description of the initial state in terms of reduced density 
matrices. 

In this paper, we will assume that the initial statistical ensemble is specifled by the 



corresponding statistical operator f?(to)- In thermal equilibrium ^(tg) is given by Eq. ( |2.5| ), 
but if the initial state deviates from thermal equilibrium, the problem of a representa- 
tive statistical ensemble becomes nontrivial. The knowledge of the Hamiltonian does not 
suffice to determine gito) and, generally speaking, one has to consider the past evolution 
of the system at times t < tp- Clearly this is not the problem which we intend to study. 
Therefore, it is reasonable to construct the initial statistical operator using a physical 
information about the system at t = tg. It is known from nonequilibrium statistical 
mechanics that a large class of nonequilibrium ensembles can be derived from Jaynes' 
maximum entropy principle [^, which is a straightforward extension of Gibbs' en- 
semble method to nonequilibrium systems. Without going into detailed discussion of this 
principle and its wide use in theory of irreversible processes (see, e.g., [ P^ pB|), we only 
recapitulate here the main ideas which will be of importance in our subsequent consider- 
ation. 

Let us assume that the state of a many-body system at t = is characterized by 
a set of observable quantities (state parameters) which are the mean values {PmY^ of 
some dynamical variables Pm- Following Jaynes' principle, the corresponding statistical 
operator, usually referred to as the relevant statistical operator, is found from a maximum 
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of the entropy S = —Tt{q' Ing'} for given (-Pm)*° — Tr|f)'Pm|, where g' is a trial 
statistical operator. /^From the extremum condition 6S = one immediately finds that 
the relevant statistical operator has the form 

^(to) =e-^/Tr{e-^} (2.6) 

with 

S = Y: KPm. (2.7) 

m 

The Lagrange multipliers are to be evaluated from the self-consistency conditions 

(Pj*« = Tr{f,(gP™} (2.8) 

for given (Pm)*". These conditions play the role of nonequilibrium equations of state. The 
operator 5* given by Eq. ( |2.7| ) may be called the entropy operator because its expectation 
value determines the entropy of the initial ensemble [^. Clearly, the equilibrium distri- 



butions are a special case of the relevant statistical operators where the Pm correspond to 
additive integrals of motion. Note also that the definition of the nonequilibrium entropy 
as the entropy in the relevant ensemble leads to a natural extension of thermodynamic 
relations to nonequilibrium states characterized by macroscopic fluxes, partial equilibrium 



in subsystems, a nonuniform distribution of local thermodynamic quantities, etc. |2^. A 
more detailed description of many-particle correlations is achieved by the use of operators 
Pm corresponding to reduced density matrices ^i2...„- Here we will not go into a discus- 
sion of the physics contained in the description of many-particle correlations in terms of 
statistical operators (|2.6|) (see, e.g. |2^, |2^). The only fact necessary for our purposes is 



that each of the relevant dynamical variables Pm in the entropy operator (p.7| ) can always 
be written as a cluster decomposition in terms of the field operators. Thus, the general 
form of the entropy operator is 

^ = Y.Ti j dri---dr'k \k{r'i . . . r^, . . . rk)ij\r'k) ■ ■ ■ ^"^(rj i){r^) ■ ■ ■ ?/'(rfc). (2.9) 
fc>i 

The quantities A^(r']^ . . . r^, . . . r^) are some functions of the Lagrange multipliers in 
Eq. ( |2.7| ). To find these functions, one has to specify the relevant dynamical variables 
Pm- In writing Eq. ( |2.9|) , we have assumed that the initial statistical operator commutes 
with a particle number operator. For superfiuids and superconductors, however, particle 
number non-conserving terms must be included 

The entropy operator (|2.9|) can be represented as a sum S = S'^ + S', where is 
bilinear in the field operators {k = 1), and S' collects all higher-order terms. In what 
follows, we shall illustrate the general theory assuming the one-particle density matrix 
^^(^o) two-particle density matrix ^'i2('^o) be independent state parameters. In 

that case the entropy operator has the form 



S = S° + S' = J dr^ dr[ X^{r[, r^) ^\r[)tlj{r^) 



+ ^ I dr^ dr^ dr[ dr'^ X2{r[r2, r^r^) ij\r2)ip\r[)t/j{r^)ij{r2), (2.10) 
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where the functions and are determined by the equations of state 
{r^gM\r[) = TT{^\r[)^p{r,)g{to)}, 
('^i^2l^i2(^o)IK'^2) = TT[ilj\r2)ip\r[)^{r^)'^{r2) g{to)} . 



[2.11] 



Generahzation to other entropy operators would be quite straightforward. 

For the relevant ensembles, it is possible to develop a Green's function formalism which 
provides a way of calculating many-particle averages and solving nonequilibrium equations 
of state by means of a diagram technique To summarize some basic relations and 
definitions (see also Appendix A), let us introduce the "evolution operator" 



(2.12) 



where x, x' are real variables, and define dynamical variables A in the "Heisenberg picture" 

Ah{x) =U{0,x)Au{x,0) =e^^ie-^^. (2.13) 

The one-particle and the n-particle thermodynamic Green's functions are then defined as 
[from here on the symbol (. . .) means averaging with the initial statistical operator ( |2.6|) ] 



^(1,1') 



(2.14) 



...n,l'---n') = i-iy (jj^il) ■ ■■^^{n)^i{n') ■ --^U^'))) , (2.15) 

where the labels (k) and {k') indicate respectively (r^,X;i.) and (r^,x'^). is the "chrono- 
logical" ordering operator with respect to the values of x; for Fermi systems, has 
the usual sign convention for permutations of the field operators. In the special case of 
thermal equilibrium, the entropy operator is S = (3 {Ji — (iN^ , so that the thermody- 
namic Green's functions ( p.l4| ) and ( p.l5| ) go over to the well-known Matsubara-Green's 
functions ||2^, if one introduces a new variable t = (3x corresponding in a sense to the 
imaginary-time evolution. 

The one-particle thermodynamic Green's function ^(1, 1') is a function of the difference 
and satisfies the Dyson equation (see Appendix A) 



fj^' dl" {G,\l, I") - /C(l, 1")} Q{l\ I') = 5(1 - 1'), 



where 



^0-^(1,1') 



d 



5{x^-x\] 



(2.16) 



(2.17) 



is the generator of "free evolution", and /C(l, 1') is the thermodynamic self-energy. The 
value of Xq may be prescribed arbitrarily. A special choice of this parameter will be 
discussed below. Equation ( p.l6|) and its adjoint can be written in the matrix form 



(2.18) 
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where the multiphcations are defined as matrix muhiphcation with respect to single- 
particle quantum numbers plus integration over x in the interval from 
The identity matrix is /(1, 1') = 5(1 — 1'). 



III. The "Mixed" Green's Functions 



In kinetic theory, the thermodynamic Green's functions are auxiliary quantities which 
can be used for evaluating characteristics of initial correlations. To study the evolution 
of the system, one has to consider real-time Green's functions. Within the standard 
^], the one-particle real-time Green's functions of interest can be put into 



formalism 18 



a matrix Green's function defined on the Keldysh contour C (see Fig. 1): 
G(l,l') = 



^++(1,1') ^+-(1,1') 
G-+(l,l') G-(l,l') 



(3.1^ 



where (1) = {r^.tT^), {!') = {r[,t[), and the field operators are taken in the usual Heisen- 
berg representation with the Hamiltonian H. The components of G(l, 1') are the causal 
Green's function g'^{l,l'), the anti-causal Green's function g°-{l,l'), and the correlation 
functions g"^{l,l'), which are defined as 



(3.2) 
(3.3) 



where TJr {T^) is the chronological (anti-chronological) time ordering operator and the 
averages are calculated with the initial statistical operator ^(tg)- The matrix Green's 
function (|3.1| ) can also be written in a compact form 



G(i,i') = -^(r^(^H(i)^l.(i'; 



(3.4) 



where is the path-ordering operator on the Keldysh contour C. The real-time Green's 
functions have the symmetry properties 



-g^'\l\ 1) 



(3.5) 



which are valid for any initial statistical operator ^(to)- 

The diagonal parts (with respect to time) of the correlation functions g are of special 
importance in kinetic theory due to relations 



^^^^(1,1') 



ti=t'^ 



(rjf?<(ti)|r'), ^<7>(l,l') 



{ri\Ql{h)\r'), 



where 



{r^\g^{ti)\r') = 6{ri - r[) + r] g^{ri,r[;t) 



(3.6) 

(3.7) 
(3.8) 
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are the nonequilibrium one-particle density matrices in the r-representation. 

The crucial point in the standard real-time Green's function formalism is that, for a 
non- correlated initial state, the function (|3.4|) obeys the Dyson equation on the contour 
C. Although in most applications the existence of the Dyson equation is not considered 
at all, this is a nontrivial fact because it implies that the Green's function G{1, 1') has the 
unique inverse function G~^{1, 1') on C. The original derivation of the Dyson equation for 
the nonequilibrium real-time Green's functions was based on a perturbation expansion in 
the interaction picture and Wick's theorem EO]. More recently it was pointed out by 



Kremp et al. that, in order to justify the existence of the Dyson equation for G{1, 1'), 
one has to specify the initial (or boundary) condition for the two-particle Green's function 

G(2)(12,1'2') = (-.)^(T4^^(l)^^(2)V^t^(2')^l,(l'))). (3.9) 

For instance, it may be the condition of weakening of initial correlations in a distant 
past [0: 



lim G'(12,1'2') t,=t, =0(1,1') 0(2,2') +r]G{l,2')G{2,l'). (3.10) 



t' =t' 



Clearly this boundary condition implies that the Keldysh contour is deformed such as tg 
— oo and the whole past evolution of the system must be involved, except when at t = tg 
the initial state of the system is non-correlated and, hence, the boundary condition ( |3.10|) 
may be replaced by the analogous initial condition of complete weakening of correlations. 
Since, in general, the evolution starts from some correlated state, the statistical operator 
g(to) does not admit Wick's decomposition, so that the existence of a Dyson equation 
for 0(1, 1') cannot be justified by a perturbative expansion in H'. Nevertheless, using 
the cluster decomposition of initial correlations in terms of the reduced density matrices 
f?i2-n('^o)' it is possible to develop a modified diagram technique for the matrix real- 
time Green's function and study some general properties of equations of motion for its 



components [15| - |Tq|. We will follow, however, the idea advocated by Danielewicz [18 



(see also Wagner's work |T9|). We introduce a matrix Green's function G on the extended 
contour C_ from Fig. 2 involving the real-time evolution and the "imaginary evolution" 
governed by the entropy operator S as in the case of the thermodynamic Green's functions. 
The structure of G_ is more complicated than that of the real-time Green's function G 
defined on the Keldysh contour. Nevertheless, by going to the "interaction picture" on 
the extended contour, it is possible to obtain for G an expression which is quite analogous 
to the expression in the standard real-time formalism with averaging over a non-correlated 
statistical ensemble (see Appendix B). Then, upon applying Wick's theorem to every term 
in the series expansion of the Green's function G in H' and S', the Dyson equation can 
be derived. In Appendix B we show that the above procedure can be performed on a 
contour C_ (Fig. 2) with an arbitrary value of the parameter Xq. We will use this fact to 
simplify the formalism by setting Xq = 0, which leads to the contour displayed in Fig. 3. 
This choice has the advantage that now G has the minimal number of components. In 
what follows the extended contour C_ will always be assumed to have the form shown in 
Fig. 3. 
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Let us start with some definitions. First of all, we introduce the variable ^ = {t, x) 
that specifies a point on the contour C_ from Fig. 3 and define integrals along the contour 
by the rule 

d^F{0= dtFit,0) - dtFit,0) + dxF{t,,x) (3.11) 

JC_ Jto on C'+ Jto on C Jq on 

The first two integrals correspond to the chronological (C^) and the anti-chronological 
(C~) branch of the Keldysh contour C, respectively. To shorten notation, we shall use the 
underlining of functions and operators defined on the contour C_. The labels {k) in such 
functions and operators mean [k) = (r^,^;.). It is convenient to introduce the function 
5{1, 2) which plays the role of the delta function on the contour C_, i.e.. 



/ di' 5(1, i')£(i') =£(!)• 

Jc 



(3.12) 



According to the integration rule ( ^.11| ), we find that 
5(1,1') = 



S{x^ - X2) S{r^ - r^) 




1,2 gC, 

1, 2 e a 



(3.13) 



otherwise, 

where S^it^ — t[) is the delta function on the Keldysh contour |T8|, p8| : 



6c{t,-t[) = l -6{t,-t[) h,t[eC-, 
[ otherwise. 

We next introduce the Heisenberg picture for operators on the contour C_: 



(3.14) 



(3.15) 



^0 — ('^OjO) being the point at the junction of the parts C and of the contour C_ 
(see Fig. 3). The operator U_{^i,^2) defined in such a way as to describe the real-time 
evolution on C and the imaginary evolution on : 



(3.16) 



where Tq is the path-ordering operator on C_, and the "effective Hamiltonian" is given by 



21(0 



H 



-iS 



(3.17) 



With the definition ( ^.15| ) of the Heisenberg picture on the contour C_, we can construct 
the corresponding one-particle and many-particle Green's functions. We shall call these 
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functions the mixed Green's functions |2^, because they involve the real-time Green's 
functions (on the Keldysh contour C) as well as the thermodynamic Green's functions 
(on the contour C^). The one-particle mixed Green's function is defined as 

G(1,1') = -^(t^ (^^(1)^1,(1'))). (3.18) 

The factor (— z) is introduced in order to have the usual definition for the real-time com- 
ponents [cf. Eg. ( ^.41 )]. Similarly, the ra-particle mixed Green's function is defined as 

(1 . . . n, 1' ■ ■ ■ n') = {Tc fe^(l) ■ ■ ■ tni^) t^M) " " " tni^')) > • (3-19) 

The difference between the mixed Green's functions and the usual real-time Green's func- 
tions is conveniently illustrated by expressing the one-particle Green's function G_ in terms 
of its components. Recalling the definition of the contour C_ (see Fig. 3), we write 

' ^(1,1') 1,1'gC, 
^>(l,l') lea^A'eC, 
g<{l,l') leCA'eC^, ^ ' ' 

[zg{i,r) 1,1'gc,, 

where G{1, V) and ^(1, 1') are the matrix real-time Green's function ( p.4|) and the ther- 
modynamic Green's function (|2.14| ), respectively. It is significant that we also have to 
consider the cross Green's functions, ^'^(1,1'), which are constructed from Heisenberg 
operators with different types of evolution: 

^>(1,1') = g>{r,x„r[t[) = -I {i^H{r,x,)i^Ur'A)) 

G<{IA') = G<{r^t^,r[x[) = -tv {^li{r[x[)^Hirih] 
These functions satisfy the obvious boundary conditions 



G(1,10 



(3.21) 
(3.22) 



^=-(1,1': 

^=-(1,1" 



f' —f 



x;=o' 



^1 — 



^6^(l,l') 



— / 

t-i — t^n 



(3.23) 
(3.24) 



which relate them to the thermodynamic and real-time Green's functions. As we shall 
see later, the cross Green's functions play the crucial role in the theory. 



IV. Dyson Equation on the Extended Contour 



As already noted, the one-particle mixed Green's function (|3.18| ) satisfies the Dyson 
equation on the extended contour C_ (Fig. 3). This equation and its adjoint can be written 
in the form 

/ rfl"Go'(l,l")G(l",l') = 5(1,1')+ / rfl"S(l,l")G(l",l'), (4.1) 
^rfl"G(l,l")Go'(l",l') =5(l,l') + ^cil"G(l,l")S(l",l'). (4.2) 
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The operator Gg^ is defined as 



where 



Go '(1,10 



Go '(1,10 


l,VeC, 


-^^o"'(l,lO 


1,1'eC,, 





otherwise, 


\ oti 2m J 


5c;(l,lO, 



(4.3) 



(4.4) 



and Qq ^(1, 1') is given by Eq. ( 2.17|) . The components of the matrix self-energy S(l, 1') 
in Eqs (|]lD and (|]2]) will be denoted by 



Ed, 10 



S(1,10 
/C<(1,10 

/c>(i,iO 

I -^IC{1,1') 



1,1' eC, 
1 e G, 1' G C,, 
1 e C,, V G C, 
1, 1' e C,, 



(4.5) 



where Zl(l, 1') is the matrix self-energy on the Keldysh contour C [|T^, 

S = 



(4.6) 



In the last line of Eq. ([4.5|) the factor (— i) is chosen from the requirement that the 
thermodynamic component of the mixed Green's function satisfy Eqs. ( p.l8| ) with the 
self-energy /C (see below) . The Dyson equations ( [4.1|) and ( [4.2|) will serve as the basis for 
our study of short-time dynamics with correlated initial states. First we wish to discuss 
some general consequences of these equations, which do not depend on the explicit form 
of the self-energy S. For this purpose, we shall consider Eqs. ( |4.1j ) and ( f4.2| ) for different 
components of the mixed Green's function G. 



A. Real-time components of the mixed Green's function 

Suppose that in Eqs. ( |4.1| ) and ( ^I^ ) the arguments 1 and 1' correspond to the Keldysh 
part G of the contour C_. Then, since the part G^ is later along C_ than G (see Fig. 3), we 
have 

+ / dl"K,<{l,l")g>{l\l'), (4.7) 
+ / dl"9'{l,l")IC'{l",l'). (4.8) 



12 



These are still matrix equations because G{1, 1') and 1') have different components 
depending on the position of the arguments 1 and 1' on the Keldysh contour C. Taking 
the (±) - components of Eqs. (|4.7| ) and (^4.8|) , and then using expressions ( |3.1|) , ( p. 141) , and 
([4.6|) together with the integration rule 



(4.9) 



one can easily obtain a system of equations for g'^^"' and g^. It is convenient, however, to 
use retarded {g^, S^) and advanced {g^, S"^) functions instead of the causal {g*^, E"^) and 
and anti-causal {g"', S") ones. The conventional definitions [|l^, are 



9'' = 9''-9^=9^- 9\ 9'' = 9''-9^ =9^- 9% 



(4.10) 
(4.11) 



Note also that 



g^l\\, 1') = ± Q\±{t, - Q\ {^>(1, 1') - g<{\, 1')} . (4.12) 

Further manipulations with Eqs. ( [4.7| ) and ( ^4.8| ) are straightforward and similar to those 
in the standard real-time Green's function formalism (see, e.g., ||2^). The only point of 
importance is that, by definition, the cross functions Q"^ and /C*^ do not depend on the 
position of their time argument on the Keldysh contour C . Taking this property into 
account, Eq. ( |4.7| ), when written in terms oi g'^ and g^^^., becomes 



+ / rfl"/C<(l,l")e?>(l",l'), 



(4.13) 
(4.14) 



The adjoint equations follow from Eq. ( |4.8| ): 

+ I rfl"6;<(l,l")/C>(l",l'), (4.15) 

^ + 11) ^'^ = ~ ^'^ + if 1") S^/^(l", 1')- (4.16) 

Again it is convenient to go over to a compact matrix notation. We define for functions 
of t (on the real-time axis) and x the multiplication as integration over the intervals 



13 



< t < oo and < x < 1, respectively, plus matrix multiplication with respect to 
single-particle quantum numbers. Then the above equations for and g^/^ take the 
form 



;^7o-'-s«)^7< = sV + /c<^^ (4.17) 

(^7o-^ - S^/^) ^7^/^ = /, (4.18) 

g^<{g,'-j:^)=g^j:^<+g<}C>, (4.19) 

g^/^ {g,' - S^M) = I, (4.20) 
where /(1, 1') = 5(1 — 1') plays the role of the identity matrix, and 



Equations ( [4. 181 ) and ( |4.2U| ) are formally identical with equations for the retarded and 



advanced Green's functions in the standard real-time formalism [|^, while Eqs. ( |4.17|) 
and ( [4.19|) , which are modified Kadanoff-Baym equations, contain explicit contributions 



from initial correlations. Note that the terms associated with initial correlations, IC^Q^ 
and Q'^K,^, enters into Eqs. ( 4.17|) and ( [4.19 ) as source terms. On the other hand, a per- 



turbative diagram expansion [|I8| leads to modified Kadanoff-Baym equations in which 
initial correlations appear as some corrections to the self-energies E^. In Appendix C we 
show that, in fact, these two descriptions are equivalent. 

B. The thermodynamic component of the mixed Green's function 

Let us turn back to the Dyson equation ( [4.1D and take the arguments 1 = (r^, x^) and 
1' = {r',x[) on the part of the contour C_ (Fig. 3). Then we arrive at the equation 

/ dl" {g^\l,l")-}Cil,l")}gil",l')=5il-l')+ f dl"}C>{l,l")Q^{l"A')- (4.22) 

Since the values of /C^(l, 1") and , 1') do not depend on whether the argument 1" cor- 
responds to the chronological (C"*") or the anti-chronological (C~) branch of the Keldysh 
contour, the last term in Eq. ( |4.22| ) is zero by virtue of the integration rule ( [4.91 ). We 
see that the thermodynamic component of G_ satisfies the Dyson equation which is not 
connected with the real-time evolution on the Keldysh contour C. The above conclusion, 
however, is almost trivial because ^(1, 1') is the thermodynamic Green's function ( |2.14| ) 



with the imaginary evolution governed by S. Nevertheless, the point to remember is that 
the existence of a closed Dyson equation for ^(1, 1') means that it has an inverse function 
^~^(1, 1') on Cx and that the component /C of the self-energy ( [4.5|) is not connected with 
the real-time evolution of the system. 
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C. The cross components of the mixed Green's function 



To complete our general discussion of the Dyson equation for G, we now derive equa- 
tions of motion for the cross components, Q^. Again we turn to Eq. ( [4.1| ) and take 1 G C, 
1' G C^. Then we obtain 

/ dl"{GQ\l,l")-^il,l")}g<il",l')=i [ rfl"/C<(l,l")e?(l",l')- (4.23) 

Here the argument t^ may be assigned to either of the two branches of the Keldysh contour 
G . It can easily be checked that in both cases we get the same equation. In the matrix 
form, it reads 

(g^^ - S«) g< = z/C< g. (4.24) 

This equation describes the evolution of ^^(1, 1') = g'^ir^t^^, r[x'^) with respect to its real- 
time argument. To derive the equation describing the "imaginary" evolution of we 
use Eq. ( [4.2| ) where we take 1 G C and 1' G G^. This gives 

/ 1") [g^\l'\ I') - /C(l", 1')} = I I 1") /C<(1", 1'). (4.25) 

Again, the time argument ti may be taken on either of the two branches of the Keldysh 
contour. Assuming for definiteness that t^ G C"*", we get the equation which is written in 
the matrix form as 

g^ (g^^ - K) = ig"" 1C\ (4.26) 



The analogous procedure can be repeated for ^^(1, 1') = g^ {r^Xi,r[t\). As a result we 
have two equations 

g>[g,^~i:^)=igic>, (4.27) 

{g^^ -K)g>= t}C>g^, (4.28) 

which describe the real-time and imaginary evolution of g^ . 

D. Relation between S and the two-particle mixed Green's function 

To obtain explicit expressions for the components of the self-energy E, one has to 
specify the form of the Hamiltonian H and the entropy operator S. If these operators are 
given by Eqs. ( |2.1| ) and (|2.10|) , then the matrix self-energy S can be expressed in terms 



of the two-particle mixed Green's function G^"^^ (see Appendix D): 

E(l, 1') = 17] I d2 dl" d2" dV" V{12, 1"2") (1"2", 1"'2+) G"^ (1'", 1') 

= 17] I d2dl"d2"dl"'G~\l,l"')G}^\l"'2-,l"2")Vil"2'\l'2), (4.29) 
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where the symbols k'^ = {rj^,^'^) and = (rj^,^^) denote points infinitesimally later 
(earlier) on the contour C_ than k = {rj^,C,k)- The "interaction matrix" V_ on the contour 
C_ is defined as 

r r(12,l'2') 1, 2, V, 2' e C, 
V{12, 1'2') = I z V(12, 1'2') 1, 2, 1', 2' G C^, (4.30) 
[ otherwise, 

where 

Vil2, 1'2') = {r,r,\V\r[r',) 5cih - h) ^cih - A) - (4.31) 
V(12, 1'2') = -\,{r,r„ r[r',) 6{x, - x,) 6{x, - x[) 6{x[ - x',). (4.32) 

Note that the matrix V_ has the symmetry property 

]/(12,l'2') =Z(21,2'1'), (4.33) 
which refiects the fact that H and S are self-adjoint operators. 

V. The Problem of Ansatz for Correlated Initial States 

A precursor of kinetic equations can be obtained from the generalized Kadanoff-Baym 
equations ( [4. 17] ) and ( [4.19| ) for Subtracting these equations, we have 



- + s< - s< + /c< - /c>. (5.1) 



The first step to a kinetic equation for the one-particle density matrix is to take the 
diagonal part of Eq. (|5.1j ) with respect to time variables. Then the left-hand side can 
exactly be expressed in terms of Qi{t), while the right-hand side still contains double-time 
correlation functions g'^. In order to convert Eq. { f).l\j into a closed equation for gi{t), 
one has to find g"^ as functionals of the one-particle density matrix. In other words, this 
is the point where the problem of ansatz arises. The simple ansatz proposed by Kadanoff 
and Baym (KB ansatz) works only in situations where the quasiparticle picture is 
adequate to describe the evolution of the system and memory effects can be neglected. In 
succeeded years efforts were directed towards refining the KB ansatz for special problems 
of kinetic theory. As already mentioned in Introduction, an important result along this 
line has been obtained by Lipavsky et al. PT|. These authors have derived exact integral 



equations which allow, in principle, to find g^ iia terms of g-^ by iteration or, at least, to 
find corrections to the KB ansatz in a self-consistent way, including memory effects and 
dynamical correlations. It should be noted, however, that the derivation of the generalized 
Kadanoff-Baym ansatz (the GKB ansatz) in Ref. |^ is based on the boundary condition 
of complete weakening of correlations in a distant past. Therefore we have to re-formulate 
the GKB ansatz for g"^ in such a way as to take initial correlations into account. 
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Another new feature of the short-time kinetics is the appearance of the cross Green's 
functions in Eq. (|5.1|) . There are two possible approaches to these functions. The first is 
to use Eqs. ([4.24|) and ( [4.271 ) as independent equations together with a kinetic equation 
for Qiit). The second is to ehminate these functions by means of some ansatz. In this 
section we will derive exact equations which allow one to formulate a generalized ansatz 
for the cross Green's functions in just the same way as for the correlation functions g'^. 

We start with the correlation functions g^ and follow the line of reasoning which is close 
to that of the work of Lipavsky et al. In order to preserve causality, it is convenient 
to introduce two auxiliary correlation functions (denoting only the time arguments) 

Glih, h) = e{h - h) g'^{h,t2), Glih, t2) = 9{t2 - h) g^^ihM). (5.2) 
Then, for ti 7^ ^2? we have 

g^'^it^M) = G%ih, h) + G\{h, h). (5.3) 

Strictly speaking, this relation does not determine g^ for ti = ^2 because G^/^ are discon- 
tinuous at this point. One can complete, however, the definition of g^ at ti = ^2 by the 
limits ti ^ ^2 =t which give the same result. Since the representation (|5.3|) will always 
be used in integrals over time, the above refinement is irrelevant. 

Let us now derive the equation of motion for G^. Differentiation with respect to ti 
gives 

' gf^ Glih, t2) = - 12) g'''{t2, t2) + d{h - ^2) ^ ^ g^^hM)- (5.4) 

The time derivative in the last term can be eliminated with the help of the Kadanoff- 
Baym equation ( |4.13| ). Using also the representation ( |5.3| ) in the time integral, Eq. ( |5.4| ) 
becomes 



^ ^ + ^ j G%{hM) - I dh S^(ti, tg) ^1(^3, t2) = t6{h - t2) g"^{t2, t2] 



+ ^^(^1-^2) j^^ dt,{T.-'{t,M) g\hM) + ^''{hM)Gl{hM)] 

+ e{u~t2)ic<{h)g>{t2). (5.5) 

A formal solution of this equation satisfying the required boundary condition at ti = ^2 
is given by 



Gl{h,t2) = ig\U, t2) g"^{t2, t2)+ r dh g^h^^) K^ih) g>{t2) 

Jt2 

+ dhj^ dUg\UM){^'^"{hM)g\UM) + ^''{hM)G\{u,t2)]. (5.6) 



This is not of course an explicit expression for G\ but a rather complicated integral 
equation. Among other things, it contains the advanced correlation functions and the 
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cross Green's function. We can derive an analogous equation for following the same 
procedure as before. Differentiating G^(ti,t2) with respect to t2 and then eliminating the 
time derivative dg'^{ti,t2)/dt2 with the help of Eq. ([4.15|) , after some algebra we obtain 

G'5(ti,t2) = -^/"(ti,ti)/(ti,t2) + r dhg<{h)lC>{h)g^{hM) 

+ i "^^^ i '^^A9''{tiM)^''"{UM) + Gl{tiM)^^{UM}9\hM)- (5.7) 

In comparison with equations for G^/^ derived by Lipavsky et al. ||2ll], the above equations 
contain new terms (the second terms on the right-hand sides) which stem from initial 
correlations. It should also be noted that the self-energies S'^ and S"^ are modified in 
the presence of initial correlations. Physically, the explicit contributions from initial 
correlations as well as the corresponding corrections to the self-energies are expected to 
go to zero in the limit to — oo, since a many-particle system "forgets" details of its 
initial state. 

Before discussing Eqs. ( ^.6|) and (|5.7] ), let us derive analogous equations for the cross 
Green's functions. First of all we note that the solution to Eqs. ( |4.24| ) and ( [4.27| ) is 



g< = ^g^}c<g + g<^^, g> = ig]c>g^ + g^,^, 

where g^om satisfy the homogeneous equations 



^hom 0) 



'horn 



9^' - 



0, 



(5.9) 



which are to be solved with the boundary conditions (|3.23| ) and (|3.24| ) for g'^{ritQ,r[x']) 
and g^iriXi^r'^to). Using the properties of the retarded and advanced Green's functions, 



lim fi(^/^(riti,r2t2) = =F^'5(n - ^2) 

ti^i2±U 



(5.10) 



and denoting only the time arguments, we may write the solution to Eqs. (|5.9| ) in the 
form 

QLmiti) = w'^ih^to) g-'ito), gLmih) = -^Q''ito) g\h,t,). (5.11) 

Inserting these results into Eqs. (p.8|) and going to a more transparent notation, we obtain 



g<{ti) = ig'^ih^to) g^ito) + 1 r dt^ g\h, ^2) /c<(t2) g. 

Jto 

g>(ti) = -ig> (to) /(to, h) + 1 1 dt^ g ic> (t2) g^it2, ti). (5.12) 

Jto 

It is interesting to note that substitution of these expressions into Eqs. ( p. 6] ) and ( |5.7| ) 
allows one to eliminate the cross Green's functions, since the ^^(to) are determined by 
the initial ensemble, as seen from the boundary conditions ( 3.23|) . However, this is not 
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a complete solution of the problem because the cross Green's functions enter into the 
self-energies /C^. 

The GKB ansatz follows from Eqs. ( ^.6|) and (|5.7|) if only the first terms on their 



right-hand sides are retained. In that case, using Eqs. ( |3.6| ) and (|573|), we obtain 

9^ih,t,) = T] [g\h,h) Qlih) - Ql{h)g^{h,t^)) . (5.13) 
Recently the same ansatz has been recovered by approximate solution of the quantum 



hierarchy for reduced density matrices |TT|. The GKB ansatz is exact in the Hartree-Fock 
approximation for the self-energies and removes some defects of the conventional KB 
ansatz |^T], |HT|. It should be particularly emphasized that the GKB ansatz is consistent 
with the exact relation (^4.12|) and, consequently, preserves the correct spectral properties 
of microscopic dynamics. Nevertheless, the range of validity of the GKB ansatz is not 
yet known with certainty because it is very difficult to make a general estimate of the 



last terms in Eqs. (|5.6| ) and ( pTZl) . Qualitative physical arguments suggest that these 
terms may be neglected if there exist two well-separated time scales characterized by the 
"collision time" and the quasiparticle lifetime r. Then the GKB ansatz provides a 
reasonable approximation for the correlation functions in the collision integral if r ^ r^, 
i.e., not far beyond the quasiparticle picture. Iteration of the integral equations (|5.6|) 



and ( |5.7| ) around the GKB ansatz allows one to find G\i^ as functional of the diagonal 



parts of the correlation functions, g^{t, t), to some order in r^/r. Note, however, that such 
a procedure can have consequences for some important properties of the kinetic equation, 
say for the conservation laws. Later we shall return to this point. 

In direct analogy to the GKB ansatz for real-time correlation functions g'^, one may 
expect that keeping only the first terms on the right-hand sides of Eqs. ( p.l2| ) is a reason- 
able ansatz for the cross Green's functions if the description of the system does not go far 
beyond the quasiparticle picture. Thus, the simple ansatz for the cross Green's functions 
reads 

g<{t)=tg''it,to)g^ito), G^{t) = -iG>{h)g\t^,t). (5.14) 

The advantages of this ansatz for the cross Green's functions are the same as those of 
the GKB ansatz ( p.l3[ ): (i) it preserves causality and describes fading memory associated 
with quasiparticle propagation, (n) effects of strong external fields can be included in 
a consistent manner, (iii) it is not connected with a quasiclassical approximation. The 
validity of this ansatz, as well as the validity of the GKB ansatz for g^ , must be checked 
in each specific situation by using the exact integral equations ( p. 61) , ( |5.7|) , and ( |5.12|) . 
These equations can also be used to find corrections to g^ and by iteration. 



VI. T-MATRix Approximation 

So far we have discussed general consequences of the Dyson equation for the mixed 
Green's function, which do not depend on the form of the Hamiltonian H and the entropy 
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operator S. In order to proceed beyond these formal results, we will consider the special 
case where the matrix self-energy S is given by Eq. ( [4.29|) . Then we can obtain different 
approximations for the components of the self-energy by determining approximate forms 
for the two-particle mixed Green's function G*^^-* . 

In this paper we restrict ourselves to a much used approximation of the many-particle 
theory, namely the T-matrix approximation. In the context of real-time dynamics, this 
approximation is quite suitable for treating short-range interactions between particles. 
Another important point is that the T-approximation satisfies general criteria for the 
conservation laws [^]. The applicability of this approximation to the thermodynamic 
component of G^^-* depends on whether the initial correlations can be described in terms 
of binary correlations. For simplicity, we shall assume that the T-approximation may be 
applied to all components of the two-particle Green's function G^^-* appearing in Eq. ( [4.29|) , 



although in some cases it would be more reasonable to use different approximations for 
its real-time and thermodynamic components. 

A. The T-matrix on the extended contour 

In the T-approximation, the two-particle Green's function is obtained by performing 



a ladder- type summation in the particle-particle channel [0, |28[|, which leads to the 
equation 

G(2)(12, 1'2') = {G(l, l')G(2, 2') +r]G{l, 2')G(2, l')} 

+ iG{l, 1") G(2,2")Z(1"2", 1"'2"') G('ni"'2"', 1'2'), (6.1) 

where repeated arguments are summed over single-particle quantum numbers and inte- 
grated over the contour G shown in Fig. 3. The T-matrix on the contour G is introduced 
through equations 

T(12,1'20 =H12,1'20 +iHl2,l''2'')G(l'',2''')G(2'',2''')T(l'''2''',l'2'), (6.2) 

T(12,1'20 =]/(12,l'20 +iT(12,l''2'0G(l^l''0G(2^2''')Z(l'''2''',l'2'). (6.3) 

In the condensed matrix notation, Eqs. ( |6.1| ) - ( |6.3| ) read 

= (G®G),^ + z(G®G)ZG(2), (6.4) 

T = V + tV{G®G)T, (6.5) 

T = V + iT{G(^G)V, (6.6) 

where ® stands for the direct product of matrices, and the index "ex" indicates the 
symmetrized [rj = 1) or anti-symmetrized {t] = —1) matrix [cf. Eq. ( |6.1| )]. 

Equation ( |6^ ) can be solved by iteration to yield the relation between the two-particle 
mixed Green's function and the T-matrix: 

G^^^ = {G^G)^^ + t{G®G)T{G®G), (6.7) 
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where the symmetrized (anti-symmetrized) 2Z-niatrix is denoted by 

T(12, 1'2') = T(12, 1'2') + r]T{12, 2'1') = T(12, 1'2') + r/T(21, 1'2'). (6.8) 

Recalhng the definition ( |4.30|) of V_, it is easy to see from Eqs. ( |6.5|) and ( |0|) that 
the components 21(12, 1'2') of the T-matrix are not zero if only the arguments 1 and 2 
correspond to the same part of the contour C_; this is also true for the arguments 1' and 
2'. Taking this into account, the T-matrix is expressed as 

T(12,1'2') 1,2 gC, 1',2'gC, 
T('19 1'9"i — </ '^^(12, 1'2') 1,2 gC, 1',2'GCa;, , s 

<i ^>^^2,1'2') 1,2 G a, 1',2'gC, ^ ' 



zr(12,l'2') l,2Ga, l',2'Ga 



X ■ 



The real-time T-matrix, T(12, 1'2'), has four components which can be put into a 2 x 2 
matrix 

T = 1^ ^— ) = ^> ) ' (6-10) 

where the superscripts (+/— ) indicate the branch of the Keldysh contour for the pair 
of time arguments. In the last line of Eq. (|6.9| ), T is the T-matrix associated with the 
"imaginary evolution" on Cx- Finally, we note that there are two cross components of the 
T-matrix, T"^, which relate the time evolution and initial correlations. 

The real-time components of Eqs. ( |6.5| ) and ( |6.6| ) yield for T the equations 

T = v + iV{G®G)T + iV{g^ ®g^)r\ (6.11) 

T = V + tT{G(^G)V + tT<{g> (^g>)V, (6.12) 

which differ from the analogous equations in the standard Green's function formalism by 
the last terms arising from initial correlations. These terms contain the cross T-matrices 

and , hence, Eqs. (|6.11|) and ( |6.12|) are not closed. Thus, though our main interest is 
with the real-time T-matrix defined on the Keldysh contour, we have to consider equations 
for the other components of the full T-matrix. 

By writing Eqs. (pD and (pj) for T(12, 1'2') with 1,2 G and 1', 2' G G^, we find 
that the "thermodynamic" T-matrix satisfies the equations 

T = v + v{g(»g)T, T = v + T{g®g)v, (6.13) 

where no real-time quantities appear, as it must be. 

In order to derive equations for the cross T-matrices, we make use of Eq. (|6.2|) , where 
1, 2 G C^. and 1', 2' G G . Then, in the matrix notation, we obtain 

= v{g ®g)T- -v{g'^ ®g'^)T. (6.14) 

Similarly, taking 1, 2 G C and 1', 2' G G^ in Eq. (U) yields 

T< = T<{g®g)v-T{g<®g<)v. (6.15) 
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The above equations can be solved formally for the cross T-matrices to give 



r> = - [I -v{g^ g)r^ v {g^- ® g^-) t, 

T< = -T {g< ® g<) v[i-{g^g) v]~' , (e.ie) 

where / is the identity matrix. To simplify these expressions, we make use of Eqs. ( |6.13|) , 
from which it follows that 

[i-v{g® g)]-^ v = v[i-ig®g) v]-' = r. (e.i?) 

Now one sees that the cross T-matrices ( |6.16| ) can be written in the form 

= -T{g^ ®g^)T, r> = -t {g> ®g>)T. (e.is) 

These expressions allow us to eliminate the cross T-matrices in Eqs. ( p.llD and ( |6.12|) . 
Then we get for the real-time T-matrix the following equations: 

T = V + iv [{G ^ G) - {g^ g')T {g> ® g>)] t, (6.19) 
T = v + iT[{G®G) - {g^ ® g^) T {g> ® g>)] v. (6.20) 

As well as being compact, they illustrate the physics of short-time processes in the presence 
of initial correlations. Let us for a moment suppose that the second terms in square 
brackets are omitted. Then Eqs. ( |6.19| ) and ( |6.20D will describe binary collisions, {G ® G) 
being a two-particle propagator in a medium. This picture of binary collisions is adequate 
at the stage of the evolution where initial correlations have been damped. The structure 
of the full Eqs. (|6.19|) and ( |6.20|) suggests that initial correlations may be regarded as an 
additional source of scattering with the corresponding T-matrix represented by T and the 
propagators represented by the cross Green's functions. It should be noted, however, that 
the "correlation scattering" described by T and the dynamical scattering described by 
the interaction matrix V proceed concurrently, which can be seen by iterating Eqs. ( |6.19|) 
and (lOol) . 

It is convenient to reformulate Eqs. ( |6.19| ) and ( |6.2CI| ) in terms of quantities which 
reflect more clearly the two-particle character of the T-matrix approximation. First of 
all, recalling the definition ( [4.31| ) of \^(12, 1'2'), we check by iteration that Eqs. ( |6.19| ) 
and (|6.2CI|) imply the following structure of T(12, 1'2'): 

T(12,1'2') = {nr2\T,2{tut[)\r[r',)6c{ti-t2)6c{t[-t',). (6.21) 

This formula may serve as a definition of an operator Tj^gl'^i) ^'i) which will be referred to 
as the two-particle real-time T-matrix. Similarly, we introduce the two-particle thermo- 
dynamic T-matrix, 7^2(^15^1)5 the interaction operator, V^2j through the relations 

T(12,1'2') = {r,r2\r,2{x,,x[)\r[r'2)6{x,-X2)6{x[-x'2), (6.22) 

V{12, 1'2') = {nr^lV.yA) 6c{h - h) S^h - t\) 6c{t[ - (6.23) 
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A comparison with Eq. ( [4.3 1| ) shows that the matrix elements of coincide with the 
interaction amphtudes in the Hamiltonian (|2.1| ). Equation ( |6.19|) is written in terms of 
the above operators (denoting only the time arguments) as 

T,,{t, t') = Vi2 Scit -t') + ij dt" F12 {G'i2(t, t") - g<,{t) T12 g>,{t")}T^2{t", t'), (6.24) 

where we have introduced the real-time and cross Green's functions describing two inde- 
pendent particles: 

Gi2(t, = ^1(^,0^2(^,0, (6.25) 

GUt) = Qi{t)<32{t)- (6.26) 



It goes without saying that Eq. (|6.24|) is a matrix equation since each of the arguments, t 



and t', can be taken on either of the branches of the Keldysh contour C . Therefore it makes 
sense to introduce two-particle operators Tj^^, T^2-> ^^"^ ^12, which are defined in the time 
interval <t < 00. Again, it is more convenient to use the retarded {T^) and advanced 
(T{|) operators instead of the chronological and anti-chronological ones. Following the 



usual rules p8[ and recalling that functions Q'^ii) do not depend on whether t is taken 
on the branch or on the branch C~ of the Keldysh contour, Eq. ( |6.24| ) is converted 
into the set of equations 

Ti^M(t, t') = S{t -t')+z r dt" Vu F,T{t, t") T,T{t", t'), (6.27) 

J to 

TUt. t') = I £ dt" dt"' r,^(t, t") {gUt", n - GUt") T12 GUn}T^2{t"'. f), (6.28) 

where the two-particle Green's functions (7^2 F^^^ are defined as 

gUt,t')=9Ht,t')gUt,t'), (6.29) 
FS^{t,t') = ±9[±it - f)] {gUt,t') - g<2{t,t')} . (6.30) 

Equation ( |6.27| ) describes effective two-particle interactions in a medium. It is similar to 
the analogous equation in the standard Green's function formalism [0, ^ and goes over 
to the latter equation in the limit to — > — cxo. Thus, the influence of initial correlations on 
the T-matrices Ti2^ manifests itself only through the short-time behavior of the propa- 
gators F^^ . On the other hand, Eq. ( |6.28| ), which can be viewed as a generalized optical 
theorem, contains an explicit contribution from initial correlations, the second term in 
braces. 

It remains to write the cross T-matrices (|6.18|) in terms of two-particle functions. 
According to the definition of the matrix multiplication, we have 

roo POO 

T<{t) = - dt'T^2it.t')GUt')^i2. ^i2(^) = -/ dt'r,2gut')T^2it'.t)- (6.31) 

■Ito Jto 
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The "thermodynamic" T-matrix, 7^2 5 is to be evaluated from Eqs. ( |6.13| ). In that case 
we have to deal only with the problem of correlations in the initial state of the system. 

B. Expressions for the self- energies 

Substitution of the two-particle mixed Green's function ( |6.7D into Eqs. ( [4.29| ) yields 
the expression for the self-energy S in terms of the T-matrix: 

1') = ir]f{12, 1'2') G(2', 2+) = ir]G{2-, 2') T(12', 1'2). (6.32) 

For its components [see Eq. ([4.5])], we have 

S(l, 1') =zr/f(12,l'2')G'(2',2+) = zr/ ^(2", 2') T(12', 1'2), (6.33) 

IC'^^l,!') = tr]T^12,l'2') {2' ,2), (6.34) 

/C(l, 1') = ri t(12, 1'2') g{2', 2+) = r/ 6?(2-, 2') f{12', 1'2). (6.35) 

The above expression for /C, combined with Eqs. ( p.l8| ) and ( |6.13| ), leads to a closed 
equation for the thermodynamic Green's function Q. However, our prime interest here is 
with formulae ( |6.33| ) and ( |6.34D which have a direct relationship to the real-time dynamics. 
Going to the components S and T.^^^ of the real-time self-energy, we infer from Eq. ( |6.33|) 
that 

J:^<{t,t')=^1^ Tr f 1^ (t, t') g>; {t',t), (6.36) 

(2) 

Sf/^(t, t') = Tr [f^l^it, t') g<it', t) + T< (t, t') g^'^t', t)} . (6.37) 

(2) 

We now turn to Eq. ( |6.34| ). The cross T- matrices can be eliminated by means of Eqs. ( |6.31| ). 
Then we get 

}Cm = -zr/ / dt' Tr {T,^,{t, t') g<,{t') T,, g>{t)} , 

"'to (2) 

lCl{t) = -ir^ / dt'jAg<{t)T,^gi,{t')T^^{t\t)}. (6.38) 

"'to (2) 

In the special case that the entropy operator contains the one-particle term only, we 
have = and, consequently, the cross components of the self-energy vanish. 

Now Eqs. (p7|) -( ^?20D , (^^ , (^l27|) , together with the expressions for the self- 
energies given in this section, provide a closed description for the evolution of the system 
in the T-matrix approximation. Initial correlations enter into play through the thermody- 
namic Green function g and the thermodynamic T-matrix 7^2? which are to be calculated 
from Eqs. ( p.l8| ) and ( |6.13| ). Although this scheme is self-consistent, it can only be re- 
alized by numerical calculations. To gain a physical picture of short-time dynamics and 
understand the role of initial correlations, it is reasonable to go over to a simpler descrip- 
tion based on a kinetic equation for the one-particle density matrix. In the next section 
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we will apply the mixed Green's function formalism to the derivation of explicit kinetic 
equations. 

VII. Kinetic Equations 

For simplicity, we assume the initial state of the system to be spatially homogeneous. 
Generalization to inhomogeneous systems is straightforward but adds of course algebraic 
complexity. In the homogeneous case, it is convenient to work in momentum represen- 
tation where all matrix elements are defined with respect to normalized single-particle 
quantum states \p) = |p, ex) in a volume fl with periodic boundary conditions. As usual, 
the limit Q oo must be taken at the end of the calculations. By going from the field 
operators to the creation and annihilation operators, aj^ and a^, we find that the one- 
particle density matrices Qf{t), given in the r-representation by Eqs. ( p.7|) and ( p.8|) , now 
have the elements 

where /^(t) = {(ipHit) ap^{t)) is the one-particle distribution function and 

m = ^+vU{t). (7.2) 
For Fermi systems {j] = —1), fp(t) is the distribution function for holes. 

A. Short-time kinetic equation in the T -matrix approximation 

Our starting point is the generalized kinetic equation (|5.1| ). Taking its diagonal 
part with respect to time variables and the p-variables, the left-hand side reduces to 
r]dfp {t)/dt. It is convenient to eliminate retarded and advanced functions on the right- 
hand side with the aid of Eq. ( |4.12| ) and the relation 

= s'^(t,t')± e[±{t-t')]{j:>{t,t')-j:^{t,t')}, (7.3) 

where E^(t,t') is a singular term in the self-energies. In the spatially homogeneous case, 
T,^ does not contribute to the diagonal part of Eq. , so that we arrive at the equation 

^ /p, (t) = V /* dt' {g;^ (t, t') s;^ (f, t) + s;^ (t, to g;^ (f, t) } 
-V f dt' {g;^ {t, t') T.;^ {f, t) + s;^ it, t') g;^ (f, t) ] 
+v {ic;^ it) q;^ it) - ^< it) ic;^ it) } , (7.4) 

where we have taken into account that the Green's functions and the self-energies are 
diagonal matrices with respect to the one-particle p-variables. 
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Equation ( |7.4| ) is still exact. We now take for the self-energies approximate expres- 
sions ( |6.36|) and ( |6.3^ ). Using also the optical theorem (|6.28|) , the kinetic equation can 
be written in the form 

i Lit) = i^f\t) + i(<;\t) 



(7.5) 



with the collision integrals 



P2 P2 



P1P2 



(7.6) 



The two-particle collision operators, /{^''(t) and Ii2'{t), are defined as 



riC), 



to 



^ dt' {g>,{t, t') {Tg gt, tS) {t\ t) + (fi^ g<, T,i) (t, t') g>,{t' , t)} 
dt' {gt,it,t') [Tggl^f^^{t\t) + {f^2 9l2T£){t.t')gUt\t)] , 



(7.7) 



/i?(t) 



+ 



i ( ^12 



to 



^ 12 ^12 



dt' { {gh - 912) it, t') {Tg T,i) {f, t) 

+ {t^2 C12 T^il) it, t') {gt2 - 9u)it', t) } 
where we have introduced the two-particle time correlation matrix 

C^2it,t') = gmT^2Quit')■ 



(7.8) 



(7.9) 



In writing Eqs. ( [7.7|) and (|7.8| ), we have followed the convention that matrix multiplication 
like ABC implies integration over all intermediate time arguments. 

The structure of /if''(t) is similar to the structure of quantum Boltzmann-like collision 
operators in the T-matrix approximation [IT^ 113, I2H]. The new ingredient is the oper- 



ator /f^(t) which involves the effects of initial correlations through the matrix Ci2it,t'). 
However, this is not the end of the story, since Eqs. ( |7.7| ) and ( [7.8] ) contain double-time 
correlation functions gi2it, t') which are to be expressed in terms of the distribution func- 
tion fp(t) for Eq. ( [7.5|) to become a closed kinetic equation. We also note that the time 
evolution of the cross correlation functions Guit) Eq. ( |7.9| ) must be specified. We shall 
make use of Eqs. (|5.13| ) and (|5.14|) to complete the kinetic equation. Then we have the 
following expressions for the two-particle quantities of interest: 



where 



9i2it,t') = 9i2it,t')QUt')Q2it') + QUt)Q2it)9Ut,t'), 
Suit) = -9i2it, to) Quito), Ghit) = -Ghito) 9t2ito, t), 

gTit,t')=g^'\t,t')g?\tX) 



(7.10) 
(7.11) 

(7.12) 
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are two-particle propagators. We show in Appendix E that Eqs. ( [7. Ill ) provide a conve- 
nient way of writing the time correlation matrix (|7.9|) : 

CuM = gUt,to)Xi2ito) 9t2M. (7.13) 

Here Xui^o) is the initial two-particle correlation matrix Xui^)^ "which is defined as the 
correlated part of the two-particle density matrix: 

Xl2W = ft2W-(ftWf?2W)ex- (7-14) 

In momentum representation, the matrix elements of gi2{t) are given by 

{PiP2\Qi2it)\p'iP'2) = (7-15) 

After inserting the expressions ( [7.10| ) and ( [7.13| ) into Eqs. ( [7.7| ) and (|7.8|) , the collision 
integrals can then be written in terms of the one-particle distribution function fp{t), the 

two-particle T-matrices T^'^{t^t')^ and the propagators gp^'^{t,t'). Thus, to obtain a 
closed kinetic description of short-time dynamics, it remains to consider equations for 
T^^^ and g^^^. Due to the properties 

'9^St,t')Y = g^y,t), {p,P,\T,^,M\p[p',r = {p'A\T,i{t',t)\p,p,), (7.16) 

we may dwell on the retarded functions only. Equation for follows directly from 
Eq. ( |6.27| ) if we use the ansatz ( [7.10| ) to express F^, given by Eq. (|6.30| ), in terms of the 
one-particle density matrix. We then obtain 

Ti^(t, t') = V,, 6{t -t')~i dt" gUt. t") {gtQ^ - f?^ f?2 T^iK^", t')- (7.17) 



Finally, the equation for g^(t,t') follows from Eq. ( 4.14 ) in momentum representation. 
Since we consider a homogeneous system, we have 



' a 



<(t,0 = 5(t-0 + 1 dt"E;^{t,t")g^^{t",t'), (7.18) 



where = |pp/2m is the free-particle energy. The self-energy is then found in 

terms of fp, T^^^, and gp^"^ by using Eq. ( |6.37D , the optical theorem ( |6.28[ ), and the 
ansatz (|5.13|) . We will not give this somewhat lengthy expression for . 

B. Short-time kinetic equation in Born approximation 

To demonstrate the principal features of the collision integrals in the kinetic equa- 
tion (|7.5| ), we shall consider the operators (|7. 7|) and ( [7.8|) in Born approximation for the 
T-matrices. For weak interaction, Eq. ( 7.17|) can be solved by iteration. The result in a 
matrix notation up to second order in reads 

T^i^(i,0 = V,,6{t-t')-tV,,g^,{t,t') {g^g^ - glg^,),,V,,. (7.19) 
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The analogous result for the advanced T-matrix is 

Tu{t,t') = V,^5it-t')+tV,^{gtQ^ - 0lQ^)t9i2{t,t')V,^. (7.20) 

Now, to find the collision integrals in Born approximation, these expressions are to be 
inserted into Eqs. (|7.7|) and (|7.8| ). Note that in all terms, except for the first term on the 
right-hand side of Eq. ( [7.8D , one can take the T- matrices to first order in the interaction 
Vi2- As already discussed, the correlation functions are eliminated in favor of the 
one-particle density matrices by means of the ansatz ( |7.1CI| ). Then, recalling Eqs. ( [7. 1| ) for 
the elements of the one-particle density matrices, a simple algebra gives for the collision 
integrals 

W = - ^ / KlpiP'St, t') {ujM> - f,JMd, ' (7-21) 

P2P'lP'2 ° 

= 2EM(ml^l2Cl2(^,^)blP2)} 



P2P'lP2 



(7.22) 

The quantities W'^^^ and W^^^ play the role of the transition probabilities or the memory 
functions and are given by 



where 



K%,P',P'S^^^') = iPiP2\Cuit,t')V^MP2) gp.p'S^'^t) {p[p'2\V^MP2)- (7.26) 
Note that the memory function ( |7.23| ) has the symmetry property 

<.l,;.^(M') = <i,,,,(M'), (7.27) 

while this is not the case for the memory function ( [7.24| ). 

A few comments should be made here about the results ( [7.21| ) and ( [7.22| ). The term 
I^^^ is a Boltzmann-like collision integral, where the memory effects are involved through 
the retarded and advanced Green's functions. Collision integrals of this type are used ex- 
tensively in the delayed quantum kinetics (see, e.g., and references therein). The extra 
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collision integral, collects explicit contributions from initial correlations and contains 
two different terms. The first term in Eq. (|7.22|) is linear in the interaction and does not 
depend on the one-particle distribution function. Lee et al. [jTO) were the first to derive 
the analogous term for a weakly coupled low-density gas. Recently the first-order correla- 
tion term in I^^^ was re-derived within the density operator theory |12| and the real-time 
Green's function formalism in which the initial correlation where specified in the form 
of a cluster expansion |2^. The remaining integral term in Eq. ( [7.22|) has a "gain- loss" 
form and gives an account of a relaxation process caused by initial correlations. To the 
authors knowledge, this term was not considered before in short-time kinetics. Physically, 
it can be relevant for describing the quasiparticle formation at the first stage of evolution. 
The most interesting thing about the "gain-loss" correlation term in Eq. ( |7.22| ) is that 
it depends on the nonequilibrium one-particle distribution functions, but its structure is 
quite different from the structure of the Boltzmann-like collision integral ( |7.21| ). Note 
that, for Fermi systems, the combination of one-particle distribution functions in this 



term, f„ f 



P2 



f f 



f f 

Jp^Jp2 



1 ~ fp^)^^ ~~ fpji nothing but the difference between 



'P2 



particle-particle and hole- hole excitations, i.e., the Pauli blocking factor. 



C. Retarded Green's function 



To complete the kinetic equation ( [7.5|) , we need to evaluate the retarded Green's 
function from Eq. ( [7.18|) . First of all we find the self-energy E^^ (t, t'), given by Eq. ( |6.37|) , 
in Born approximation for the T-matrices. To eliminate the correlation function g"^, we 



make use of the GKB ansatz ( |5.13|) . Then we obtain 



Kit.t') = s?f (t) S{t - t') + Y.KpS^^t')9iit',t), 

P2 



(7.28) 



where 



S?f(t)=E(ml^l2blP2)/p,W 

P2 



is the Hartree-Fock term and we have defined 

1 I ~ 2 

\p, it^t') = -J2\ {PlP2 1 \P'1P2) I 9p[p'^ (t, t') [fp^ fp>^ fp,^ - fp^ /p. /p. ) ^, 



2 

P1P2 



+ V {PlP2\^l2 Cuit^t') 1^12 biPs)- 

Insertion of Eq. (|7.28| ) into Eq. (|7.18|) leads to the integro-differential equation 



where 



Ep^it)=ep^ + J:f^{t) 



(7.29) 



(7.30) 



I -i^p.w) 9^,M = S{t-t')+^j'^dt''Ap^p^{tXMSt,t'r9;^S^ (7.31) 



(7.32) 
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is the re-normalized particle energy. Equation ( |7.31D should to be solved together with 
the kinetic equation for the one-particle distribution function, but it seems to be a very 
difficult task, even using numerical methods. We therefore will construct an approximate 
solution to Eq. ( |7.31| ), valid for a weak interaction. 

First we note that the retarded Green's function t') can always be written in the 
form 

g;^^{t,t') = -te{t - t') exp - <fp^(t,t')} , (7.33) 

where 

= /*^t"^p,(n = ^p,(t-o+/'^t"s?f(n, (7.34) 

and $p (t, t') is some complex- valued function. It satisfies the obvious initial condition 
t) = and has a simple physical interpretation. Its real part determines the quasi- 
particle damping, while the imaginary part can be associated with higher-order corrections 
to the quasiparticle energy. Inserting the expression (|7.33|) into Eq. (|7.31|) , one obtains 
for $p^(t,t') an integro-differential equation, which can then be recast into the following 
integral equation: 



P2 



where 



(7.35) 
(7.36) 



^p,P2 (^, = ^p, (^ > t') + (i, • 

Since the right-hand side of Eq. ( |7.35| ) is already of second order in the interaction, the 
leading approximation for $p(t,t') can be obtained by setting $ = in the integrand. In 
the weak coupling case, the imaginary part of $p(t,t') is small compared with Up{t,t'), 
so that our interest is in the real part of $p(t, t') describing the quasiparticle damping. 
With Eq. (17:301) , we find 



{PlP2\yi2\p'lP2) 



r,_((,(') = Re {*,_((,(')} 



- E 

2 V, 

P2P1P2 



X cos 



E / dr^ j " dr^ Re { 

P2 



{PiP2\yi2Ci2ir^,r2)V^2\PiP2)} , (7.37) 



where 



^'^PlP2'P'lP2 



(7.38) 



The first term in the expression ( |7.37|) is long-lived because a substantial contribution to 
the inner integral is made by the region (r^ — Tg) of the order of the collision duration 
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time and, consequently, the first term grows with {t — t') if (t — f ) ^ r^. The behavior 
of the second term depends critically on the time rj„ required for damping of the initial 
correlations. Since the matrix elements of Ci2{t,t') according to Eq. (|7. 13|) are oscillating 
functions of t — and t' — t^, the second term in Eq. ( [7. 371 ) is expected to have a plateau 
at (t — to) ^ and (t' — tg) ^ r^„. The interplay between the "quasiparticle" and 
"correlation" contributions to Tp{t, t') at the initial stage of the evolution can be of interest 
in short-time kinetics. However, a systematic analysis of this point requires more detailed 
information on the interaction operator and the initial two-particle correlation matrix 

Xl2(io)- 

To summarize, we may conclude that a reasonable approximation for the Green's 
functions in Eqs. ( [7.13| ), ( [7.25| ), and ( [r.26| ) is given by 



-_-^e{t- f) exp { ~iuj^^ {t,t')- (t, t') } . 
This retarded Green's function satisfies the equation 

- E^, (t) + {t, t')) gl^it, t') =6{t- t') 



(7.39) 



(7.40) 



with the time-dependent damping 

d_ 

di 



E 



2 ~T , Jt 

P2P1P2 



dr 



(PiP2l'^i2b'iP2) cos[Aa;p^p^ (t, r)J [fpjp'jp'^ - fpjp'jp'^ 



~vJ2f, dr Re {e^-^1^2 {p^p^ I V,2C,2 (t , r ) \p^P2 ) } • 
P2 *' 

Equation (|7.40|) can be considered as an approximate version of Eq. (|7.31|) . 

D. Memory effects and energy conservation 



(7.41) 



The expressions ( [7.39| ) for the retarded and advanced Green's functions with u^^ (t, if:') 
and ^p_^it, t') given respectively by Eq. (|7.34|) and Eq. ( [7.37|) complete the kinetic equation 
since the memory functions in the collision integral, as well as the matrix Ci2{t,t'), can 
now be written in terms of the one-particle distribution function and the initial correlation 
matrix Xui^o)- The most important feature of the short-time kinetic equation is the non- 
Markovian structure of the memory functions. Let us write, for instance, the explicit 
expression for the memory function (|7.23|) using Eq. ( [7.39| ). We have 



{PlP2\Vl2\PlP2) 



(t,t') 



COS 



Au;, 



P1P2 



where the following designation is used: 



^v.vM = {t. t') + r^^ (t, t') + r^, (t, t') + r^, (t , t') 



(7.42) 



(7.43) 
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The structure of the memory function ( [7.421 ) is typical for the memory functions used in the 
delayed quantum kinetics (see, e.g., 0), but the damping function Tp(t,t') is commonly 
taken in a simplified form rp(t,t') = 7p(t — t') with a constant 7p. In some sense this 
corresponds to the asymptotic form of our result ( |7.37| ). It should be noted, however, that 
at the first stage of evolution such an approximation seems to be very rough; we see from 
Eq. ( |7.37D that the behavior of T^it, t') is rather complicated and involves the contribution 
from initial correlations. Another source of memory effects is the oscillating cosine term in 
Eq. ( [7.42D , which replaces the energy-conserving delta function in Boltzmann-like kinetic 
equations. We see that to take into account the full structure of the memory functions, 
even in performing the numerical evaluation of short-time quantum kinetic equations, 
is very complicated. Nevertheless, there are some necessary conditions which must be 
checked for any kinetic equation. The most important issue is whether or not the kinetic 
equation obeys the conservation laws. The problem of conservation laws in the delayed 
quantum kinetics is not trivial because of strong memory effects and the influence of initial 
correlations. We shall discuss the conservation laws for the kinetic equation, where the 
collision integrals are given by Eqs. ( [7.21|) and (|7.22|) . Since there are no problems with 
the conservation of the total number of particles and the total momentum, even for the 
general form of and g^, we restrict our discussion to the conservation of energy, which 
is a serious problem in kinetic theory. 

Since we are working in the Heisenberg picture, the average energy of the system is 
written as £{t) = {H{t)). In our case the Hamiltonian has the form ( p.l|) . Going over 
to momentum representation and taking into account that the system is assumed to be 
spatially homogeneous, we find that 

= ES,4W + lj:iPlP2\yi2012it)M, (7.44) 

Pi P1P2 

where gi2{t) is the nonequilibrium two-particle density matrix with the elements given by 
Eq. ( |7.15|) . It is convenient to separate the correlated and non-correlated parts of gi2(t) 
by using Eq. ( [7. 14|) . Then the average energy (|7.44|) takes the form 

The quantity S^^j^{t) is the average kinetic energy; this is just the first term in Eq. ( [7.44| ). 
The Hartree-Fock contribution to the average energy, S-^p{t), is given by 

^Hf(^) = ^ E {PlP2\Vl2\PlP2) fpS^)fp,it). (7.46) 

P1P2 

Finally, the term 

ScoiM = ^ E {PiP2\V^2Xi2it)\PiP2) (7.47) 
P1P2 

may be called the correlation energy. Energy conservation implies that the time derivative 
of S{t) is zero. Having a kinetic equation for fp^it), one can calculate 

|(^kin(^)+^^HF(^)) =E^P,(^)^%^ = E^P,(^)^p,(^), (7-48) 

Pi Pi 
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where Ep^(t) is the re-normahzed particle energy, Eq. ( |7.32| ), and Ip_^it) is the collision 
integral. Note that the above relation is exact. Thus, the necessary condition that a 
kinetic equation is consistent with the energy conservation is that the sum on the right- 
hand side of Eq. ( [7.48| ) can be represented as the time derivative. If this condition is 
fulfilled, then we have 



EE, 

Pi 



Pi 



dt 



(7.49) 



This formula allows one to evaluate the approximate correlation energy from the kinetic 
equation. 

It is interesting to check the condition ( |7.49| ) for the short-time kinetic equation, in 
which the collision integral Ip^{t) is the sum of the collision integrals I^^\t) and I^\t) 
given by Eqs. ( |7.21| ) and ( [7.22| ). The crucial point is the choice of the retarded and 



advanced Green's functions entering into these collision integrals. It turns out that the 
necessary condition for the energy conservation is satisfied if gp^{t,t') and g^^{t' ,t) are 
taken in the form ( [7.39|) , where the damping function Vp^{t,t') is set to be equal to zero. 
The corresponding manipulations are straightforward but somewhat lengthy, so that we 
only explain the main points. Multiplying Eqs. (^) and ([^221) by (t) and then 
summing over p^, the products of the type Ep(t)gp(t,t') or g^{t\t)Ep{t) appear in each 
term. These products are eliminated with the aid of Eq. (|7.40| ) and the corresponding 
equation for gp(t',t) (recall that Tp(t,t') = 0). It can then be shown that the 5-functions 
do not contribute to the sum J2p Ep(t)Ip(t). The remaining terms form the time derivatives 
of the contributions to the correlation energy. In this way, we arrive at Eq. (|7J|) with 
the correlation energy 



^corr (^) 



E 



dt' 



{PiP2\Vi2\p'iP2) sin[A^ y , (t, f) 



Plpip'lp'i ° 



^ (/pi /p2 •^P'l •^P2 fpJp2^P'l^P'2)t' 
+ \T.^^{{PlP2\yi2Cl2it.t)\PlP2)} 



P1P2 



E [dt'H<Up'.M] {fp'Jp'. - fp'Jp'X' ■ 



(7.50) 



PlP2P'lP2 ° 



Here the retarded and advanced Green's functions entering into Ci2i't,t) and X^^\t,t') 
are understood to be taken in the approximation described above (with Tp{t,t') = 0). 
An attractive feature of the expression (|7.5CI| ) is that it gives the exact result for the 
initial correlation energy. To show this, we recall Eq. ( |7.13| ), from which it follows 
that C^2i^O'^o) ~ Xui^o)- Thus we see that, indeed, the initial values i^corr('^o) given 
by Eqs. ( p'.47| ) and ( |7.50| ) are identical. 



Some previously obtained results follow from Eq. (|7.50|) as special cases. If all the 



contributions from initial correlations are omitted and one sets ui {t,t') = e Jt — t'), the 
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expression ( |7.5UD reduces to the correlation energy calculated from the non-Markovian 
Boltzmann equation [^. The second term in the correlation energy (linear in the in- 



teraction) agrees with the density operator result |T2|, Q obtained by truncating the 
quantum hierarchy for the reduced density matrices. The last contribution to the corre- 
lation energy (|7.50|) stems from the "gain-loss" term in the collision integral (|7.22 ). 



It may appear tempting to improve the expressions for the collision integrals and 
the correlation energy by taking the retarded and advanced Green's functions in the 
form ( [7.39| ) with {t,t') 7^ 0. This means that one goes over to the memory functions 
involving effects of the quasiparticle damping [see, e.g., Eq. (|7.42| )]. Such a procedure 
is often used (with an empirical exponential damping factor) in numerical calculations 
based on non-Markovian Boltzmann kinetic equations and in many cases is crucial for 
the stability of the results. Unfortunately, the inclusion of the quasiparticle damping into 
the memory functions leads to collision integrals which do not conserve the total energy 
of the system. The origin of this difficulty is easy to see. We recall that the above proof 
of the energy conservation rests heavily on Eq. ( |7.4CI| ) for the retarded Green's function. 
If %^{t,t') = 0, then the product Ep^{t)g^^{t,t') can exactly be expressed in terms of 
the time derivative of the retarded Green's function. For a finite damping, however, this 
is not the case, and the right-hand side of Eq. (|7.49|) involves extra terms violating the 
energy conservation. Clearly, the contribution of these extra terms to the energy balance 
equation is determined by the time behavior of the retarded Green's function. Strictly 
speaking, the total energy is conserved only if the product 'yp^(t,t')gp^{t,t') is equal to 
zero for all t — t'. At the initial short-time stage of evolution 7p^(t,t') is completely neg- 
ligible, so that the total energy is conserved. The real difficulties appear in the long-time 
asymptotic behavior of the collision integrals, where 7p^(t, t') is close to its stationary 
value. Then the decay of g^ {t, t') must be sufficiently fast for the energy non- conserving 
terms in the collision integral to be negligible. For the electron scattering with optical 
phonons, Haug and Banyai |Q proposed a model retarded Green's function which van- 
ishes faster than the slowly decreasing function with the constant quasiparticle damping. 
Their numerical results show that the collision integral with the improved gp{t,t') con- 
serves the total energy much better in comparison to the non-Markovian collision integral 
with the exponential decay of the memory function. It is more natural, however, to de- 
scribe the transition from the short-time dynamics to the asymptotic Boltzmann regime 
by using approximate self-consistent solutions of the equation for the retarded Green's 
function. For instance, the result (|7.41| ) shows that the effective quasiparticle damping has 
a non-Markovian structure. Furthermore, it depends on the nonequilibrium one-particle 
distribution functions and involves contributions from initial correlations. Nevertheless, 
even if one uses the self-consistent expression for gp{t,t'), the problem of energy conser- 
vation for all time intervals still remains open since any quasiparticle damping function 
7p(t, t') (except for ■jp{t, t') = 0) leads to the energy non-conserving collision integral. This 
is not surprising, because the extra terms in the energy balance equation are formally be- 
yond the Born approximation used in the derivation of the kinetic equation. Thus, to 
obtain the collision integrals which involve the quasiparticle damping and are consistent 
with the energy conservation, one has to calculate the T-matrices in the full collision 
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operators ( |7.7| ) and ( |7.8| ) keeping terms of higher-order in the interaction. 

Concluding this section, we would like to touch upon another important aspect of the 
short-time kinetics. Suppose that the evolution of the system starts from the equilibrium 
state. Then a kinetic equation must lead to the obvious result dfp{t)/dt for all t > t^. In 
other words, the collision integral must be equal to zero in thermal equilibrium. It is well 
known that the stationary solution of the Markovian Boltzmann equation corresponds 
to the equilibrium Fermi or Bose distribution function for quasiparticles. If, however, 
the collision integral takes account of initial correlation and memory effects, the question 
about the stationary solution of the kinetic equation becomes nontrivial. We have shown 
that the kinetic equation with the collision integrals ( 7.21|) and (|7.22|) conserves the total 
energy (if Tp{t, t') = in the memory functions) and that the initial value of the correlation 
energy given by Eq. ( [7.50| ) coincides with the exact result. Due to this fact, it is reasonable 
to expect the kinetic equation to have the correct stationary solution. Nevertheless, we 
should show strictly that the collision integral vanishes in thermal equilibrium, which is 
a stronger condition than the conservation of energy. In Appendix F we calculate the 
collision integrals ( [7.2 Ij ) and (|7.22|) up to terms of second order in the interaction using 
the equilibrium two-particle correlation matrix x^i2^ ■ The resulting kinetic equation has 
the form 



|4W = - E \{piP2\vMp'2) 



dt' 



COS 



Aw, 



PlP2'PiP2 



+ E \{PlP2\Vl2\p'lP'2) 



P2P1P2 



2 sm 



PlP2'PlP2 



^PlP2'P'lP'2^^-^^ "^^^^ 



^P[ - ^P'2 



(7.51) 



where is the particle energy in thermal equilibrium, including the Hartree-Fock con- 
tribution, and the function is defined as 



'^PlP2'P'lP'2^^'^^^ fpifp2'^P'l'^P'2 fpxfp2^P'Jl: 



P2 



(7.52) 



The first term on the right-hand side of Eq. ( [7.51| ) is the non-Markovian Boltzmann 
collision integral ( |7.21| ), where the quasiparticle damping is omitted. The second term 
is due to the equilibrium correlations. In this term fjf^^ is the equilibrium one-particle 
distribution function. Since, in equilibrium. 



^^P,P2Ap'2(t^t') = {e^^ + E^^ - E^, - E^,) {t - t' 



(7.53) 



it can easily be seen that = /^'^'^^ is a stationary solution of Eq. ( |7.51| ). Physically, this 
means that in equilibrium the changes in the one-particle distribution function caused 
by collisions and correlations exactly cancel each other. Lee et al. [|lOl were the first 
to demonstrate this fact for a weakly interacting low-density classical gas. The kinetic 
equation ( [7.51| ) shows that the same situation holds in quantum systems. It should also 
be noted that in the latter case the Hartree-Fock corrections to the particle energies have 
to be included, in order to assure the conservation of energy. 
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The fact that in equihbrium the colhsion and correlation effects cancel each other 
is analogous to the fluctuation-dissipation theorem and suggests that correlations and 
fluctuations are closely related. The interplay between collisions and correlations becomes 
particularly clear if one neglects the Hartree-Fock terms in the memory functions. Then 



both terms on the right-hand side of Eq. ( 7.51 ) can be combined to give 



\{PlP2\yi2\p'lP2) 



P2P1P2 



dt' cos 



'Pi 



'P2 



X 



(7.54) 



Neglecting here the correlation term we recover the well-known Levinson equation |35 
The advantage of Eq. (|7.54| ) in comparison with the Levinson equation is that it has 
the right equilibrium solution. Although the correlation term was introduced in the 
equilibrium form, Eq. (|7.54|) can be considered as a rather simple and a reasonable kinetic 
equation describing short-time quantum kinetics in the presence of initial correlations. 



VIII. Conclusion 

We now summarize the main features of the outlined approach to short-time quantum 
kinetics: 

1) The formalism is applicable to a wide range of correlated initial states described by 
nonequilibrium statistical thermodynamics. 

2) Initial correlations and the time evolution are incorporated on the same footing 
through the many-component ("mixed") Green's functions. The essential point is that 
the introduction of the mixed Green's functions makes only minor changes in the standard 
real-time Green's function method and allows the use of the Dyson equation for correlated 
initial states. The existence of the Dyson equation for the mixed Green's function leads 
to some exact relations which play a crucial role in the derivation of kinetic equations. 
In particular, the generalized Kadanoff-Baym ansatz can be formulated for short-time 
processes in the presence of initial correlations. 

3) The well-known approximations, originally developed within the standard real-time 
Green's method, say the T-matrix approximation, are conveniently extended to short-time 
kinetics with correlated initial states. 

4) The method given in this paper is non-perturbative in external fields due to the 
self-consistency of the generalized ansatz for the real-time and cross correlation functions. 
This is important when considering a direct effect of a strong external field on collisions. 
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APPENDIX A: 

Thermodynamic Green's Functions 



Here we briefly discuss some properties of the thermodynamic Green's functions in 
terms of which initial correlations are described. Owing to the structure of the evolution 
operator ( p.l2| ), these properties are in many ways similar to those of the Matsubara- 
Green's functions In particular, the one-particle thermodynamic Green's function 

Q{1, 1') = Q{riXi,r[x'i) depends on Xi and x'l through the difference follows 
at once from the cyclic invariance of the trace in Eq. (|2.14|) . Note also that ^(1,1') is 
defined for — 1 < ( ) < 1, which is a necessary condition for convergence of the 

trace. Clearly this condition is satisfied if the "imaginary" Heisenberg picture ( |2.13| ) is 
introduced for x^ < x < Xq + 1, where Xq is an arbitrary parameter. Then it is easy to 
verify that ^(1,1') has the property 



(A.i) 



which is a generalization of the well-known Kubo-Martin-Schwinger boundary condition 
for the equilibrium Green's functions ^ . 

Perturbation expansions of the thermodynamic Green's functions can be constructed 
by introducing the interaction picture 



Aj{x) 



-,xS' 



Ae 



(A.2) 



where Xq < x < Xq + 1, and 5**^ is the entropy operator for a non-correlated state [cf. 
Eq. (p.lO|) ]. The interaction picture evolution operator is defined as 



h{j{x, x') = 
Then, using the identity 



^xSO^-{x-x')S^-x'S° 



exp - / dx"S'j{x"] 



W,(0,Xo)e-^°W,(xo + l,0), 



(A.3) 



(A.4) 



which follows directly from the definition of Uj{x,x'), the one-particle thermodynamic 
Green's function ( p.l4|) can be written in the form 



T^(l(,ix, + l,x,)Ml)4il'] 



(A.5) 



{Uj{xo + 1,Xo))q 

where the symbol (. . .)o stands for averages calculated with the statistical operator 

,°(to) = e-^7Tr{e-^°}. (A.6) 

Expression (|A.5|) has the structure typical for a diagram technique. Since 5*° is bilinear in 
the field operators, the statistical operator p° admits Wick's decomposition and, conse- 
quently, expectation values of products of the interaction-picture field operators factorize 
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into one-particle thermodynamic Green's functions. The corresponding Feynman rules 
depend on the specific form of the "correlation" term S' in the entropy operator. 

In the diagram language, the thermodynamic self-energy JC is introduced as an irre- 
ducible part of Q. It can also be defined without reference to diagrams. To do this, let 
us consider equations of motion for the field operators in the Heisenberg picture ( p.l3|) . 
They read 



d 

dxi 

d 
dxi 



til'Ai(l,l')V^H(l')+[^H(^i)>V^H(l 
cil>i,(l')Ai(l',l) + fe(xi),^l,(l) 



where dl' = dx\ dr[, and the function Xi{l, 1') is defined as 



(A.7) 



(A. 



Then, differentiating the Green's function ( |2.14| ) with respect to Xi and x[, and using the 
above equations of motion for the field operators, we arrive at the Dyson equations ( 2.18|) 
with the self-energy 



/C(l,l') 



dl" {T:{[s'^ix,),i;M] 

dl" g-\i, 1") {t: {mi") [^1,(1'), s'Ax[)] } 



(A.9) 



These relations can be used to express the self-energy in terms of higher-order ther- 
modynamic Green's functions. If the entropy operator is given by Eq. ( p. 101) , we have 
(integration over repeated arguments is implied) 



/C(l,l') = r]Vil2,l"2")g^^\l"2",l"'2+)g-\l"',l'), 
= vS-\l, 1"') g^^\l"'2-, 1"2") V(l"2", 1'2), 

where the amplitude V is given by Eq. ( ^.321 ). 



(A.IO) 



APPENDIX B: 

Interaction Picture on the Extended Contour 

Let the Heisenberg picture on the contour C_ of Fig. 2 be defined as 

AHio=m^o,OAm^,^o). (B.i) 

where the variable ^ = {t,x) specifies a point on C_ and £[(^1,^2) evolution opera- 

tor ( pj.l6| ). For the contour C_ shown in Fig. 2, the "effective Hamiltonian" 2i.{0 coincides 
with H on the Keldysh part C and with —iS on the parts C'^ and C". 
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The definition ( [B.l| ) of tlie Heisenberg picture on C_ allows us to introduce the path- 
ordered ("mixed") Green's functions, as given by Eqs. ( p.l8|) and ( |3.19|) . Due to the 



fact that the evolution over the Keldysh contour C is described by the identity operator 
U.{^o,^o) = 1, the mixed Green's functions coincide, up to a factor, with the thermo- 
dynamic Green's functions discussed in Sec. II and Appendix A, if all the arguments 
correspond to the parts C'^ and C'^ of C_. On the other hand, if all the arguments of the 
mixed Green's function correspond to the Keldysh contour, then this function coincides 
with the real-time (path-ordered) Green's function. The latter property follows directly 
from the cyclic invariance of the trace in Eqs. (|3.18|) and ( p.l9|) . 



To formulate perturbation theory for the mixed Green's functions defined on the con- 
tour C_ of Fig. 2, let us introduce the "interaction" picture through the relation 

MO = U\^o,OAU\^,^^). (B.2) 

The evolution operator is 

U'{^^,Q=Tcexpl-z f'n\Od^], (B.3) 



-.0 , 



where the unperturbed "effective Hamiltonian" Ti. (0 coincides with on the Keldysh 
contour C and with —iS'^ on the parts C'^ and C'^ of C_. Then an operator in the Heisenberg 
picture (|B.1D can be written as 



where 

is the interaction picture evolution operator. This operator is represented by a path- 
ordered exponent 

1^/(^1,^2) =%exp|-2^^'?i;(0rfe| (B.6) 
with the effective "interaction Hamiltonian" 

HjiO = { ^ B.7 

Here Hi{t) is an operator in the real-time interaction picture and Sj{x) is the correlation 
part of the entropy operator in the thermodynamic interaction picture introduced in 
Appendix A. The representation ( p.6|) can be derived by solving the equation 

1/7(^1,^2) =21/(^1) £7(^1,^2), (B.8) 

which follows from Eq. (|B.5|) . 
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Let us now consider the product of Heisenberg operators AiHi^i) ' ' ' AkuHk)-: where 
the points are arranged in a certain order on the contour C_ of Fig. 2. Making use of 
Eq. ( p.4| ), we can write this product in terms of the interaction picture operators: 



(B.9) 



As a next step to perturbation expansions of the mixed Green's functions, we use the 
relation 

e-^ = f/,(eo-,eje-^°I^/(eend,eo-), (B.IO) 
which is nothing but the identity ( |A.4D since 

HAio^^-J = W/(0, X,), UriCn^, ^o) = W/(^o + 1,0), (B.ll) 

as is evident from Fig. 2. Now, combining Eqs. (|B.9|) and ( p.lO|) , the one-particle mixed 
Green's function ( ^.181 ) takes the form 





exp 








1). 


( 


^T^exp 




c_ 





:b.i2) 



where averages are calculated over the non-correlated initial ensemble described by the 
statistical operator ( |A.6|) . 

Note that the representation ( [B.12|) is valid for any value of the parameter Xq in the 



interval — 1 < Xq < 0. It is convenient, however, to take Xq 



-1 or Xq = 0. To illustrate 



this point, let us take Xq = (see Fig. 3), as in the main body of the paper. Then, 

expanding Eq. ( |B.12| ) in terms of 21/(0 ^^"^ applyi^g Wick's decomposition to each term 
in this expansion, one will only obtain the cross Green's functions ( p.21|) and ( p.22| ), in 
which the x-arguments are always later on C_ than the t-arguments. In particular, the 
expansion of the real-time component of G_ can be derived from formula 



G(l, 1') = -I 



W,(1,0)T^ exp 



J dtH'.it) ^h(I)^Ui') 



c 



(W,(1,0)), 



(B.13) 



where all the operators with the imaginary evolution are arranged to the left of operators 
with the real-time evolution. For a more general contour C_ shown in Fig. 2, one has to 
introduce not only the cross Green's functions ( p.21| ) and ( |3.22| ), but also two functions 



^>(l,l')=^>(riti,r;x;) 
^<(1,10 =^<(riXi,r^t;) 



ir] (^j/(r'/i)V^i^(riXi) 



(B.14) 
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which appear when the x-argument is taken on the part C'^ of C_. It should be emphasized, 
however, that of the four cross functions, and JF^, only two are independent of each 
other due to relations 

[g>{i, 1')]* = --F>(i', 1), [g<{i, 1')]* = --F>(r, i), (b.is) 

where (k) = (r^,-^^). 



APPENDIX C: 

An Alternative Form of Generalized Kadanoff-Baym 

Equations 



We shall discuss the connection between our equations ( [4.17| ), ( [4.19| ), and those de- 
rived by the diagram technique . First of all we note that Eqs. (|5.8| ) yield the following 
expressions for the correlation terms appearing in the generalized Kadanoff-Baym equa- 
tions: 



horn 



ic<g> = iic<gic>g^ + ic<gi^^, 

We now make use of Eqs. (|5.11|) to write 

where we have introduced the singular self-energies 

^""{tuh) = i5{U - to - 0) g^{t^) JC^ih). 
^ci^iM) = -ilC<{t,)g>{to)6{t2-to-0). 



(C.l) 

(C.2) 

(C.3) 
(C.4) 



Combining Eqs. (|C.1|) with ( p.2|) , we see that Eqs. (|4.17| ) and ([4.19|) can be rewritten in 
the form 



9o' 



S +^c)9 



where 



(C.5) 
(C.6) 

(C.7) 



are the re-normalized real-time components of the self-energy. 

Equations ( |C5| ) and ( |C.6| ) are similar to the equations derived by Danielewicz from 



diagram expansions of Green's functions |3^. In the diagram language, S are irreducible 
parts of the real-time Green's function, which begin and end with an interaction amplitude 
in the Hamiltonian H |jl8|. On the other hand, the singular self-energy T,^ C^c) 
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represented by diagrams which begin (end) with a correlation matrix and end (begin) 
with the interaction amphtude V. In our approach, such a structure of the above self- 
energies follows directly from equations of motion for Q^, Eqs. ( [4.24|) and (|4.27|) , since 
the function /C< (/C'') contains the interaction amplitude as a left-side (right-side) 
multiplier. 



APPENDIX D: 
Self-Energy on the Extended Contour 



If the entropy operator (|2.9| ) contains a finite number of the correlation terms, one can 
derive a hierarchy of equations for the mixed Green's functions, which is analogous to the 



Martin-Schwinger hierarchy in the standard real-time Green's function formalism |28 



Truncation of this new hierarchy at some order allows one to formulate reasonable ap- 
proximations for the self-energy on the extended contour C_, as much as truncation of 
the Martin-Schwinger hierarchy is used for evaluation of the self-energy on the Keldysh 
contour C. Here we will discuss briefly the connection between the self-energy E and 
higher-order mixed Green's functions in the model where the entropy operator includes 
only the two-particle correlation term, Eq. (p.lO|) . 

The hierarchy of equations for the mixed Green's functions follows directly from equa- 
tions of motion for the field operators ^ ^{1) and 'ip''^{l) in the Heisenberg picture on the 
extended contour C: 



(D.l) 



Recalling the definition (|3.17| ) of the "effective Hamiltonian" on the contour C_, the above 
equations can be transformed to (with integration over repeated arguments) 



G-,\l, 1') = H12, 1'2') f ^(2) l^{2') 1^(1'), 

t^^il') G-o\l\ 1) = Z(1'2', 12) f^(l') I^J2') ±^{2\ 



(D.2) 



where the operator G g ^ and the "interaction amplitude" V_ on the contour C_ are defined by 
Eq. ( [4.3|) and Eq. (|4.3CI|) , respectively. Now Gq^ operating on the Green's function ( |3.18|) 
from the left and right and then using Eqs. ( p.2|) results in the following equations: 



GfH, 1") G(l", 1') = 5(1, 1') + zr/ V:(12, l"2") G^^\l"2\ 1'2+), 
G(l, 1") Go'(l", 1') = 5(1, 1') + ir]G^^\l2-, l"2") V{l"2" , l'2), 



(D.3) 



These are the first equations of the hierarchy for the mixed Green's functions on the 
contour G, which is the analogue of the Martin-Schwinger hierarchy in the standard real- 
time Green's function formalism p8[. Assuming that G_ has an inverse on the contour G, 
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Eqs. ( |D.3| ) can be written as Dyson equations ( |4.1j ) and ( [4 .21 ) with the matrix self-energy 
E given by Eq. (^1291) . 

We emphasize that it has been possible to express the self-energy S in terms of the two- 
particle mixed Green's function G^^-* due to the structure of the entropy operator ( p.lO| ). 
For a more general form of the entropy operator, the self-energy E will depend on higher- 
order mixed Green's functions. 



APPENDIX E: 

The Two-P article Time Correlation Matrix 



To calculate the time correlation matrix ( |7.9| ), we use the ansatz ( 7.11 ) for the cross 
Green's functions and find that 

Let us now turn to Eqs. ( |6.13| ) for the thermodynamic T-matrix. Recalling Eq. ( |6.22|) 
and the boundary conditions ( |3.23| ), we write 



6^r2(io)^i2e^r2(io) = f'dxdx'g,^iO,x)T,^{x,x')g,^{x',0), (E.2) 
Jo 

where we have introduced the notation 

gi2{^, x') = Qi{x, x') g^ix, x'). (E.3) 

We now wish to show that the right-hand side of Eq. ( [E.2|) can be expressed in terms of 
the two-particle thermodynamic Green's function g^'^^ [see Eq. ( p. 151 )]. To that end we 
use Eq. (|6.7|) which, when written for the thermodynamic component of G^^-* , reads 

g('\i2,i'2') = (g{i,i')g{2,2')) 

\ /ox 

+ g{i, 1") g{2, 2") t(i"2", i"'2"') g{i"', i') g{2"', 2'). (e.4) 

It is convenient here to go to the matrix notation with respect to the single-particle 
quantum numbers. Defining the matrix ^^2'' (^12^25 ^'1^2) t>y 

e;(2)(i2,l'2') = {r,r,\gi'^{x,x,,x[x',)\r[r',), (E.5) 

we have 

g[2 {^1X2, X1X2) = X]^) ^2(^2? -^2))^^ 

+ dx" dx'" g^X^.x") g2{x2, X") f^2ix", X'") g^{x'", X[) g2{x"', X'2). (E.6) 
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A comparison with Eq. ( |E.2| ) shows that 

^12(^^0) ^12 Suito) = lim {Quixx, x'x) - {Q^[x, x) Q^[x, a;'))^^} • (E.7) 



Using the obvious relation 

^r^r^\ Q^yl[xx,x'x')\r'^r'^ = e{x - x') (?/'^(rix)^/'^(r2x)V'H(r2x')^]^ (r^a;')) 

+ 9{x' - x) {^Ur'^x')^l{r[x')ijHirix)ijHir2x)), (E.8) 

it is easy to verify that the order of the hmits in Eq. (|E.7] ) is of no significance, as it 
should be. In both cases we obtain the same result 



where Xui^o) is the initial two-particle correlation matrix (|7.14| ). Relations ( |E.1|) and ( [E.9| 
complete the derivation of Eq. ( [7.13| ) . 



APPENDIX F: 

Derivation of the Kinetic Equation with Equilibrium 

Correlations 

The first (Boltzmann) term on the right-hand side of Eq. ( |7.51| ) follows immediately 
from Eq. ( |7.21|) if we take the memory function W^^^ in the form ( [7.42| ) with = 0. 



Thus it remains to evaluate the correlation term given by Eq. (|7.22|) . Since in the case of 
equilibrium correlations the matrix Ci2{t,t') is proportional to and we wish to obtain 
J^*"^ to order Vj^^, we can restrict our discussion to the first term in Eq. ( |7.22|) . We start 
by evaluating the equilibrium correlation matrix 

(P1P2IX1? \PlP2) = (4^4iSiS2)eq - fil'^^fpl''^ {(^p,p[^p^p'., + VSpy./p^p[) , (F.l) 



where the symbol (■ ■ ■)^^ stands for averages calculated with the statistical operator (|2.5|), 
and f^^'^^ = {a^pap)eq is the equilibrium one-particle distribution function. The second- 
quantized Hamiltonian in momentum representation is given by 

H = H' + H' = Y.^pA.^p. + 1 ^ {PiP2\Vu\p'iP2) (^-2) 

Pl PlP2P'lP'2 

Recall that we need to evaluate the correlation matrix to first order in the interaction. 
This can be done by elementary methods by expanding the statistical operator ( |2.5|) in 
H'. It is convenient, however, to redefine the unperturbed Hamiltonian by taking the 
Hartree-Fock term into the particle energies. Then, instead of Eq. ( p.5| ) we now have 



= ^-m'+H') /Tre-^(^°+^'), (F.3) 
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where we have defined 



Pi 



PlP2P'lP2 ^1 



Pi%r 



The re-normahzed particle energies are given by 

E — ^ I ^^HF 



Pi ^P, + ^P^ = Si + II (^1^2 I '^12 \PlP2) fp^ 



(eq) 



P2 



Up to terms hnear in Ti' Eq. ( [b'.3D reads 



-eq 



-eqi 



where the statistical operator 



-eq 



Tr e 



(F.4) 
(F.5) 

(F.6) 

(F.7) 
(F.8) 



describes the ideal quantum gas with the quasiparticle energies ( [F.6|) . The symbol (■ ■ ■)q 
means the average with g'^^. Making use of Eq. ( [F.7|) , it is easy to see that to first order 
in the interaction the one-particle distribution can be replaced by the Fermi or Bose 
distribution 



(eq) 



(F.9) 



We now can evaluate the average value in Eq. ([F.l|) with the aid of Eq. ( [F.7| ). Using the 
relations 



we find that 

(PiP2lXi?|p'iP2) 



E. 



P1P2 



^p'iP'2 



1 1 (al,^al,a^^apn 



(F.IO) 



(F.ll) 



where E„ 



.p^p^ ^p^ Ep^. In the above expression the superscript (c) shows that only 
the connected part of the average must be taken. This means that, in applying Wick's 
theorem, all the creation and annihilation operators in Ti' must be contracted with the 
"free" operators. Further manipulations are very simple. After some algebra we find 



{PlP2\Xl2^\p'lP2) 



^p^P2,p[p'X^f^'"^^}) 



E. 



P1P2 



^P'rP'2 



{PlP2\yi2\p'lP'2) , 



(F.12) 



where the function JF({/}) is defined by Eq. ( |7.52| ). In writing Eq. ( |F.12| ), the temperature 
dependent factors have been eliminated by means of Eq. (|F.9|) . 

Now we substitute the expression ( p^'.12| ) into Eq. ( [r.l3| ) and take the retarded and 
advanced Green's functions in the form (|7.39|) with F^ = 0. This gives us the time- 
correlation function Ci2{t,t) which is to be used for evaluating the first term in the 
collision integral ( [7.22| ), i.e., the correlation contribution into the kinetic equation ( [7.51|) . 



45 



References 

[1] R. Zwanzig, J. Chem. Phys. 33 (1960), 1338. 

[2] I. Prigogine and P. Resibois, Physzca 27 (1961), 629. 

[3] R. Balescu, Physica 38 (1968), 98. 

[4] J. Shan, "Ultrafast Spectroscopy of Semiconductors and Semiconductor Microstruc- 
tures". Springer Series in Solid-State Sci., Vol. 115, Springer, Berlin/Heidelberg, 
1996. 

[5] H. Haug and A.-P. Jauho, "Quantum Kinetics in Transport and Optics of Semicon- 
ductors", Springer Series in Solid-State Sci., Vol. 123, Springer, Berlin/Heidelberg, 
1996. 

[6] P. Danielewicz, Ann. Phys. (N.Y.) 55 (1984), 305. 

[7] C. Greiner, K. Wagner, and P. G. Reinhard, Phys. Rev. C 49 (1994), 1693. 
[8] H.S. Kohler, Phys. Rev. C 51 (1995), 3232. 
[9] H.S. Kohler, Phys. Rev. E 53 (1996), 3145. 
[10] D. Lee, S. Fujita, and F. Wu, Phys. Rev. A 2 (1970), 854. 

[11] D. Kremp, M. Bonitz, W.D. Kraeft, and M. Schlanges, Ann. Phys. (N.Y.) 258 
(1997), 320. 

[12] M. Bonitz, "Quantum Kinetic Theory", Teubner, Stuttgart/Leipzig, 1998. 

[13] L.P. Kadanoff and G. Baym, "Quantum Statistical Mechanics", Benjamin, N.Y., 
1962. 

[14] L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47 (1964), 1515 [Sov. Phys.- JETP 20 (1965), 
1018]. 

[15] A.J. Hall, J. Phys. A 8 (1975), 214. 
[16] A.J. Hall, Physica A 80 (1975), 369. 

[17] Yu.A. Kukharenko and S.G. Tikhodeev, Zh. Eksp. Teor. Fiz. 83 (1982), 1444 
[Sov. Phys.- JETP 56 (1982), 831]. 

[18] P. Danielewicz, Ann. Phys. (N.Y.) 55 (1984), 239. 

[19] M. Wagner, Phys. Rev. B 44 (1991), 6104. 

[20] D. Semkat, D. Kremp, and M. Bonitz, Phys. Rev. E 59 (1999), 1557. 



46 



[21] P. Lipavsky, V. Spicka, and B. Veliky, Phys. Rev. B 34 (1986), 6933. 
[22] E.T. Jaynes, Phys. Rev. 106 (1957), 620; 108 (1957), 171. 

[23] E.T. Jaynes, "Information Theory and Statistical Mechanics". In: "Statistical 
Physics. Brandeis Lectures", Vol.3, Benjamin, N.Y., 1963, p. 181. 

[24] D. N. Zubarev, V. G. Morozov, and G. Ropke, "Statistical Mechanics of Nonequi- 
librium Processes", Vol. 1, Akademie Verlag, Berlin, 1996. 

[25] D. N. Zubarev, V. G. Morozov, and G. Ropke, "Statistical Mechanics of Nonequi- 
librium Processes", Vol. 2, Akademie Verlag, Berlin, 1997. 

[26] G. Ropke, Ann. Physik (Leipzig) 3 (1994), 145. 

[27] A.L. Fetter and J.D. Walecka, "Quantum Theory of Many-Particle Systems", 
McGraw-Hill, N.Y., 1971. 

[28] W. Botermans and R. Malfliet, Phys. Rep. 198 (1990), 115. 

[29] R.A. Craig, J. Math. Phys. 9 (1968), 605. 

[30] D. Kremp, M. Schlanges, and Th. Bornath, J. Stat. Phys. 41 (1985), 661. 
[31] K. Morawetz and G. Ropke, Phys. Rev. E 51 (1995), 4246. 
[32] K. Morawetz, Phys. Lett. A 199 (1995), 241. 

[33] M. Bonitz, D. Kremp, D.C. Scott, R. Binder, W.D. Kraeft, and H.S. Kohler, Journal 
of Physics: Condensed Matter 8 (1996) 6057. 

[34] H. Haug and L. Banyai, Solid St. Commun. 100 (1996) 303. 

[35] LB. Levinson, Zh. Eksp. Teor. Fiz. 57 (1969) [ Sov. Phys. - JETP 30 (1970) 362 ]. 



[36] It should be noted that, in Danielewicz's paper |T^, there is no clear distinction 



between the two different self-energies, S^andE . In Eq. (5.22) of ReL [l|], the S< 
have the same meaning as in the present paper. On the other hand, in Eqs. (H.13) 
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FIGURE CAPTIONS 



Figure 1: 

The Keldysh contour C with the lower (chronological) branch C"*" and the upper (anti- 
chronological) branch C~ 

Figure 2: 

The extended contour C_ with the real-time evolution on the Keldysh contour C and the 
"imaginary" evolution on the parts and C'^. The parameter xq satisfies — 1 < < 0. 

Figure 3: 

A special case {xq = 0) of the contour shown in Fig. 2 



